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Foreword 


This  report  summarizes  the  progress  of  four  technical  investigations  conducted  during 
CY  88.  Although  reasonable  efforts  have  been  made  to  ensure  the  reliabihty  of  the 
data  presented,  it  must  be  emphasized  that  this  is  an  interim  progress  report  and  that 
further  experimentation  and  analysis  may  be  performed  before  the  conclusions  from 
any  of  these  investigations  are  formally  published.  It  is  therefore  possible  that  some 
of  the  observations  presented  in  this  report  will  be  modified,  expanded,  or  clarified 
by  our  subsequent  research. 
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Research  for  Electric  Energy  Systems  - An 

Annual  Report 

R.  J.  Van  Brunt,  Editor 


Abstract 

This  report  documents  the  technical  progress  in  the  four  investigations  which 
maJte  up  the  project  “Support  of  Research  Projects  for  Electrical  Energy  Sys- 
tems”, Department  of  Energy  Task  Order  Number  137,  funded  by  the  U.S. 
Department  of  Energy  and  performed  by  the  Electrosystems  Division  of  the 
National  Institute  of  Standards  and  Technology  (NIST).  The  first  investigation 
covers  the  measurement  of  ions  in  the  vicinity  of  dc  high-voltage  transmission 
lines  and  in  biological  exposure  facilities  designed  to  simulate  the  environment 
around  transmission  lines.  The  measurements  of  ion  mobilities  at  atmospheric 
pressure  using  two  different  types  of  drift  tubes  are  discussed.  Results  are  also 
presented  from  an  investigation  into  the  uncertainties  associated  with  mea- 
surement of  monopolar  space-charge  densities  in  the  atmosphere  around  high- 
voltage  transmission  lines  where  dc  electric  fields  are  present.  The  second  area 
of  investigation  is  concerned  with  the  behavior  of  gaseous  dielectrics  like  SFe 
used  in  high-voltage  power  systems.  Results  from  two  significant  technical  ac- 
tivities in  this  area  are  summarized  here.  The  first  is  concerned  with  mezisure- 
ment  of  collisional  electron  detachment  probabilities  for  negative  ions  moving 
in  gaseous  SFe  and  the  implications  of  the  results  of  these  detachment  mea- 
surements on  achieving  a better  understanding  of  ion  transport  and  positive 
discharge  initiation  probabilities  in  SFe.  The  second  is  about  a new  method 
for  measuring  the  stochastic  behavior  of  pulsating  partial-discharge  phenom- 
ena. Although  this  method  is  applied  here  to  an  investigation  of  Trichel-pulse 
corona  in  gas  mixtures  like  air  and  SFe/02,  it  is  argued  that  it  has  a more  gen- 
eral application  to  measurement  of  partial  discharges  and  offers  the  possibility 
of  revealing  much  more  information  about  the  physical  bases  for  these  processes 
than  can  be  achieved  by  presently  used  measurement  techniques.  The  third  area 
of  investigation  has  to  do  with  observation  of  prebreaJcdown  phenomena  in  di- 
electric liquids.  Results  are  presented  here  from  optical  observations  of  the 
influence  of  hydrostatic  pressure  on  prebreakdown  partial  discharge  develop- 
ment and  from  measurement  of  nanosecond  impulse  breakdown  at  liquid-solid 
interfaces.  An  image-preserving  optical  delay  system  is  described  for  use  in  ob- 
taining photographs  of  improbable  fast  breakdown  events  in  liquids.  The  fourth 
area  of  research  discussed  in  this  report  is  concerned  with  electrical  measure- 
ment of  fast  transient  phenomena.  Results  are  presented  from  an  investigation 


vi 


into  the  interactions  between  two  dividers  used  simultaneously  to  measure  fast 
impulse  voltages.  A discussion  is  given  about  the  difficulties  encountered  in 
attempting  direct  comparisons  of  the  responses  from  two  impulse  dividers  used 
simultaneously.  Techniques  for  minimizing  the  effects  of  the  interactions  on 
measurement  accuracy  are  considered. 
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1.  ELECTRIC  FIELD  MEASUREMENTS 


Task  01 

Martin  Misakian,  William  Anderson,  and  Ronald  McKnight 
Electricity  Division 

National  Institute  of  Standards  and  Technology 


1.1  Introduction 


The  objectives  of  this  project  are  to  develop  methods  to  evaluate  and  calibrate  instru- 
ments which  are  used,  or  are  being  developed,  to  characterize  the  electrical  parameters 
in  the  vicinity  of  power  lines  and  in  laboratory  apparatus  designed  to  simulate  the 
power  line  environment.  Electrical  measurement  support  is  also  provided  for  DOE- 
funded  efforts  to  determine  if  there  are  biological  effects  due  to  ac  fields  and  dc  fields 
with  ions. 

The  electrical  parameters  of  interest  include  the  electric  field  strength,  magnetic  flux 
density,  ion  current  and  charge  densities,  ion  mobility  and  ion  species.  Recent  lab- 
oratory studies  have  focused  on  the  measurement  of  ion  mobilities  at  atmospheric 
pressure  using  drift-tubes.  During  the  past  year  the  evaluation  of  two  drift  tubes 
fabricated  from  insulating  cylinders  was  completed.  Excerpts  from  this  study  are 
presented  in  the  following  sections.  A manuscript  which  describes  the  full  study  has 
been  accepted  for  pubhcation  [1].  The  influence  of  a dc  electric  field  on  the  mea- 
surement of  monopolar  charge  densities  using  an  aspirator-type  ion  counter  and  the 
measurement  of  net  space  charge  density  using  a Faraday  cage  or  filter  was  also  ex- 
amined. Optimum  configurations  which  minimized  the  effect  of  the  electric  field  were 
identified  for  each  type  of  instrument.  Some  results  of  this  study  are  briefly  described 
in  this  report.  A full  account  of  the  work  will  be  published  at  a later  date  [2].  Other 
NIST  activities  in  1988  included:  three  revisions  of  a new  IEEE  draft  standard  which 
will  provide  guidance  for  the  measurement  of  dc  power  line  fields  and  ion-related 
parameters;  a site  visit  to  the  Battelle  Pacific  Northwest  Laboratories  to  characterize 
electrical  parameters  in  several  bioeffects  exposure  systems;  and  invited  conference 
talks  at  the  IEEE  Instrumentation  and  Technology  Conference  and  the  NSF  spon- 
sored U.S.- Japan  Seminar  on  Electromagnetic  Interference  in  Highly  Advanced  Social 
Systems.  Both  talks  considered  measurement  procedures  for  characterizing  electrical 
pcirameters  near  ac  and  dc  power  lines. 
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1.2  Ion  Mobility  and  Drift  Tubes 

The  mobility  of  an  ion  as  it  moves  through  a gas  under  the  influence  of  a weak  electric 
field  is  defined  as  the  ratio  of  the  ion  drift  velocity  to  the  electric  field  strength. 
Corrections  for  pressure  and  temperature  are  normally  made  to  obtain  a “reduced 
mobility.”  As  has  been  noted  previously,  the  motivation  for  measuring  ion  mobilities, 
within  the  context  of  possible  bioeffects  due  to  dc  power  lines,  is  that  mobility  values 
provide,  in  some  cases,  a rough  measure  of  ion  mass  [3,  4].  Ideally,  the  ion  species  is  of 
interest  when  considering  the  possibility  of  bioeffects  due  to  exposure  to  ions  from  dc 
power  lines,  but  sampling  ions  at  atmospheric  pressure  with  a mass  spectrometer  is 
difficult  and  requires  a complex  apparatus  [5].  Ion  mobility  spectra  can  be  measured 
at  atmospheric  pressure  with  less  complicated  instrumentation  and  ascertaining  that 
the  mobility  spectra  remain  unchanged  during  a bioeffects  study  may  be  as  important 
as  knowing  that  the  electric  field  strength  or  ion  current  density  remains  constant. 

Ion  mobility  can  be  determined  by  measuring  the  time  required  for  an  ion  to  travel 
a given  distance  in  a known  dc  electric  field.  A device  that  that  has  been  used 
extensively  for  this  measurement  is  the  drift  tube.  The  two  drift  tubes  that  were 
evaluated  are  shown  in  Figme  1 and  have  been  described  in  detail  previously  [6].  The 
construction  of  the  drift  tubes  differs  from  those  of  conventional  design  in  that  metal 
guard  rings  with  insulating  spacers  and  alignment  rods  have  been  replaced  with  an 
insulating  tube  on  which  guard  rings  have  been  coated. 

Figure  1(a)  shows  a drift  tube  designed  to  operate  in  a pulsed  time-of-flight  (TOF) 
mode.  Briefly,  ions  arriving  from  the  right  are  prevented  entry  into  the  drift  tube  by 
the  shutter  which  is  biased  shut  with  a small  voltage  supplied  to  the  outside  mesh. 
The  ions  are  produced  in  untreated  room  air  by  wires  in  corona,  using  an  apparatus 
that  has  been  previously  described  [7].  Ions  enter  the  drift  tube  and  travel  to  the 
collector  after  the  shutter  is  pulsed  open  for  a fraction  of  a millisecond.  Because  more 
than  one  species  of  ions  is  present  in  the  approximately  disk-shaped  volume  of  ions 
that  passes  through  the  shutter,  the  ions  will  separate  into  groups  with  the  same 
drift  velocity  or  mobility  as  they  traverse  the  length  of  the  drift  tube.  A shield  grid 
before  the  collector  prevents  the  induction  of  current  in  the  collector  circuit  until  the 
ions  pass  through  the  shield  grid.  The  current  induced  in  the  collector  circuit  as  the 
ions  travel  between  the  shield  and  collector  is  monitored  and  exhibits  maxima  and 
minima  as  a function  of  time  because  of  the  different  time-of-flights  of  the  different 
ion  species.  The  length  of  the  drift  region  is  variable  so  that  the  distance  the  ions 
travel  can  be  made  nominally  4 cm  or  8 cm.  By  obtaining  TOF  spectra  at  two  drift 
distances,  it  is  possible  to  eliminate  so  called  “end  effects”  when  computing  the  TOF 
of  a given  ion  peak. 

Figure  1(b)  shows  a schematic  view  of  an  ac  TOF  drift  tube  similar  in  principle  to 
that  employed  by  Tyndall  and  Powell  [8].  The  shutters  at  each  end  of  the  drift  tube 
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Figure  1.  Schematic  views  of  drift  tubes  made  of  insulating  cylinders  with  conductive 
guard  rings  on  the  inside  walls,  (a)  pulsed  time-of-flight  drift  tube  with  movable 
collector  electrode.  Metal  tabs  at  L equal  to  4 cm  and  8 cm  (movable)  provide  appro- 
priate potential  to  shield  grid,  (b)  ac  time-of-flight  drift  tube.  The  inside  diameter 
of  each  tube  is  5.65  cm. 
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are  biased  closed  and  are  opened  periodically  with  sinusoidal  voltages  applied  to  the 
inner  mesh  of  each  shutter.  Ions  that  have  TOF’s  between  the  shutters  equal  to  half 
the  period  of  the  ac  voltage  will  pass  through  both  shutters  to  yield  a steady-state 
current  which  is  measured  with  an  electrometer.  Ideally,  measurement  of  ion  current 
as  a function  of  frequency  results  in  current  maxima  at  specific  frequencies  for  ions 
with  different  mobilities.  The  mobility  of  an  ion,  K,  is  given  by  the  expression 

K = 2L^f/V,  (1) 

where  / is  the  frequency  of  the  current  maximum,  L is  the  length  of  the  drift  space 
and  V is  the  voltage  across  the  drift  tube. 


1.3  Experimental  Results 

At  the  close  of  1988,  it  had  been  determined  that  (1)  for  the  pulse  TOF  drift  tube, 
there  existed  a discrepancy  between  the  shape  and  peak  position  of  measured  TOF 
spectra  and  that  determined  from  theory  and  (2)  the  use  of  sinusoidal  voltage  wave- 
forms to  gate  the  shutters  of  the  ac  TOF  drift  tube  failed  to  resolve  ions  with  different 
mobility  in  the  TOF  spectrum  [6].  The  resolution  of  these  problems  is  described  be- 
low. 

(1)  Matching  Theoretical  and  Experimental  Spectra:  Pulse  TOF  Drift  Tube 

Figure  2 shows  fits  of  the  fastest  TOF  feature  of  complex  positive  ion  spectra  obtained 
on  different  days  for  a drift  space  of  8 cm.  The  fits  are  made  using  trial  values  of  drift 
velocities  in  a theory  described  in  detail  previously  [6].  This  theory  takes  into  account 
such  factors  as  diffusion,  shutter  pulse  width,  and  transit  time  between  the  shield 
grid  and  collector.  The  choice  of  drift  velocity  is  made  by  matching  the  calculated 
initial  rise  of  the  ion  current  with  the  experimental  one.  Ideally,  the  fit  should  be 
to  the  shape  of  the  spectrum,  but  simultaneous  matching  of  the  onset  and  peak  ion 
currents  in  figure  2 was  not  possible.  The  data  in  figures  2(a)  and  2(b)  are  obtained 
with  feedback  resistors  in  the  amplifier  circuit  [R  in  fig.  1(a)]  of  2000  and  200  MD, 
respectively. 

It  is  found  that  the  fitting  process  yields  mobility  values  that  are  in  good  agreement 
with  values  determined  by  obtaining  TOF  spectra  for  two  different  drift  distances. 
This  latter  approach  is  expected  to  yield  the  most  accmate  mobility  values  because  it 
is  free  of  end  effects.  While  there  is  agreement  between  the  mobility  values  determined 
experimentally  (1.97  ± 0.02  x 10“^  m^/Vs),  the  shape  of  the  calculated  spectrum  as 
noted  above  does  not  match  the  data  very  well  in  figure  2(a).  The  agreement  between 
mobility  values  suggests  that  the  peak  of  the  measured  spectrum  may  be  shifted  to  a 
longer  TOF  due  to  (i)  the  possible  overlapping  with  a slower  adjacent  feature,  (ii)  an 
instrumental  effect,  or  (iii)  some  combination  of  (i)  and  (ii). 
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Figure  2.  TOF  spectra  of  positive  ions  with  the  largest  mobility,  obtained  with  the 
pulse-TOF  drift  tube  {L  = 8 cm)  on  different  days.  The  data,  represented  by  the  solid 
curves,  are  fitted  with  theoretical  curves  (dashed  curves)  and  assmning  the  indicated 
mobility  values,  (a)  measurement  performed  with  2000  megohm  feedback  resistor  in 
current-to-voltage  amplifier  circuit,  (b)  measurement  performed  with  200  megohm 
feedback  resistor  on  different  day.  While  the  spectra  overlap  with  slower  adjacent 
features,  the  discrepancies  between  theoretical  and  measured  peak  TOF  values  is  due 
mainly  to  excessively  long  circuit  time  constants. 
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While  there  is  some  overlap  of  the  fastest  TOF  feature  with  an  adjacent  slower  one, 
mechanism  (ii)  appears  to  be  the  predominant  reason  for  the  discrepancy  in  peak 
position.  Because  of  the  2000  megohm  feedback  resistor,  stray  capacitances  of  a few 
tenths-of-picofarads  in  the  circuit  can  lead  to  time  constants  of  several  tenths-of- 
milliseconds  and  cause  the  observed  shift  in  the  TOF  of  the  peak  ion  current.  This 
explanation  is  supported  by  the  results  in  figure  2(b)  which  shows  improved  agreement 
between  the  calculated  fit  and  data  when  the  feedback  resistor  is  reduced  by  a factor 
of  ten.  The  calculated  and  experimental  mobility  values  are  1.92  ± 0.02  x 10“^  m^/Vs 
and  1.92  ± 0.02  x 10“‘‘m^/Vs,  respectively,  and  the  difference  in  peak  positions  is 
reduced  to  about  0.1ms.  A risetime  of  about  40/zsec  in  an  optical  link  circuit,  not 
shown  in  figure  1,  also  contributes  to  the  discrepancies  in  TOF’s  in  figure  2.  Further 
evidence  for  the  time  constant  explanation  is  provided  below  when  new  data  obtained 
using  the  ac  TOF  technique  are  discussed. 

(2)  Resolving  Ions  With  Different  Mobilities  Using  AC  TOF  Drift  Tube 

As  has  been  previously  noted  [6],  a theoretical  expression  for  a TOF  spectrum  is 
not  readily  developed  when  the  shutters  of  the  ac  drift  tube  are  gated  open  with  a 
sinusoidal  voltage.  However,  a more  serious  difficulty  is  that  at  atmospheric  pressure, 
the  frequency  of  the  gating  voltage  can  be  below  50  Hz.  At  such  low  frequencies 
and  for  the  ac  voltage  amplitudes  required  to  produce  adequate  signal  strength,  the 
duration  of  the  open  condition  can  be  several  milliseconds.  The  long  duration  results 
in  a very  broad  spectrum  making  it  difficult  to  identify  the  frequency  in  Eq.  1 . Because 
more  than  one  ion  species  is  normally  present  in  the  spectrum,  the  long  opening  of 
the  shutters  also  results  in  the  loss  of  resolution  of  the  individual  current  maxima. 

During  1988,  the  sinusoidal  waveform  was  replaced  with  a pulse  waveform  of  alter- 
nating polarity.  This  modification  permitted  the  derivation  of  an  expression  which 
describes  the  spectrum  shape  as  a function  of  frequency  of  a single  ion  species  and 
it  also  became  possible  to  resolve  ions  of  different  mobility  in  the  experimental  spec- 
tra. [1].  With  the  alternating  pulse  waveform,  the  expression  for  ion  mobility  that 
corresponds  to  Eq  1 becomes  [1] 


= FL’^jV.  (2) 

F is  the  inverse  of  the  ion  TOF  and  is  given  by 

f = [(l//p) -«]-'>  (3) 

where  fp  is  the  experimental  frequency  corresponding  to  the  peak  ion  current  and  8 
is  the  delay  in  ion  TOF  as  the  ion  passes  through  the  first  shutter  (8  is  determined 
experimentally). 

To  gate  the  shutters  of  the  ac  drift  tube  in  the  ac  pulse  mode,  the  transformers  in 
figure  1(b)  were  replaced  with  phototransistor  optoisolators  as  shown  in  figure  3. 
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Figure  3.  Schematic  view  of  circuit  which  converts  pulses  from  the  function  generator 
to  pulses  of  alternating  polarity.  Nearly  identical  pulse  voltages  are  apphed  to  the 
inner  mesh  of  each  shutter  to  maintain  the  same  electric  field  in  the  drift  space. 
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Figure  4.  Positive  ion  current  as  a function  of  shutter  frequency.  The  current  measure- 
ments are  made  with  the  ac  TOF  drift  tube,  using  alternating-polarity  pulse-voltages 
(pulsewidth  = 0.73  ms)  to  gate  the  shutters.  The  dashed  curve  indicates  the  general 
features  of  the  spectrum. 


A train  of  rectangular  pulses  from  a waveform  generator  is  converted  into  pulses 
of  alternating  polarity  by  means  of  TTL  flip-flop  NAND  gates.  The  same  voltage 
waveform  is  apphed  to  the  inner  mesh  of  each  shutter  as  for  the  sinusoidal  waveform. 

An  example  of  spectra  that  is  obtained  when  the  sinusoidal  waveform  is  replaced  by 
a pulse  waveform  of  alternating  polarity  is  shown  in  figure  4.  The  data  represent 
positive  ion  collector  current  for  two  sweeps  of  the  frequency  from  20  Hz  to  114  Hz 
in  2-Hz  increments.  Each  point  is  an  average  of  five  current  measurements  using  an 
electrometer-voltmeter  combination.  While  the  structure  in  figure  4 is  reproducible, 
in  general,  considerable  variation  of  the  peak  intensities  occurs  above  80  Hz.  The 
data  in  figure  5 represent  four  sweeps  of  the  pulse  frequency  from  85  Hz  to  113  Hz  in 
1-Hz  increments  and  is  compared  with  a calculated  spectrum  (sofld  curve)  determined 
from  theory. 

For  the  spectrum  in  figure  5,  T = 6.13  ± 0.03  cm,  fp  = 101  ± 1 Hz,  V = 2000  ± 2 V, 
and  6 = 0.27  ± 0.02  ms.  From  Eqs  (2)  and  (3),  the  value  of  mobility  is  1.95  ± 0.03  x 
10“^m^/Vs.  The  mobility  determined  from  the  theoretical  fitting  process  is  also 
1.95  ± 0.02  X lO--*  mVVs. 


1.3  Experimental  Results 
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Figure  5.  Positive  ion  current  as  a function  of  shutter  frequency  for  the  same  condi- 
tions in  figure  4.  The  solid  curve  is  a theoretical  fit  of  the  fastest  feature.  The  data 
in  this  figure  and  in  figure  2(a),  using  the  pulse-TOF  drift  tube,  were  obtained  nearly 
simultaneously. 
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The  results  in  figure  5 should  be  compared  with  the  spectrum  shown  in  figure  2(a), 
which  was  measured  at  the  same  time  using  the  pulse-TOF  drift  tube.  All  of  the 
mobility  values  agree  within  the  noted  uncertainties.  Also,  it  can  be  shown  that  the 
theories  for  the  pulse  and  ac-pulse  spectra  are  the  same  except  for  the  technique  of 
detecting  the  ion  current  [1].  The  ac-pulse  measurement  technique  involves  the  mea- 
surement of  an  average  steady-state  ion  current  at  each  frequency  or  TOF  and  circuit 
time  constant  problems  are  avoided.  Thus,  the  good  agreement  between  the  shape  of 
the  data  and  the  calculated  spectrum  in  figure  5 provides  further  evidence  that  the 
discrepancy  between  the  calculated  and  experimental  peak  TOF’s  in  figure  2(a)  is 
due  to  stray  capacitance  effects  in  the  current  detector  of  the  pulse-TOF  drift  tube. 


1.4  Conclusions 

This  study  examined  the  operation  of  two  types  of  drift  tubes  that  are  suitable  for 
measuring  ion  mobilities  at  atmospheric  pressure  in  the  presence  of  large  dc  electric 
fields.  The  fabrication  of  the  drift  tubes  from  insulating  cylinders  with  guard  rings 
coated  on  the  inside  wall  differs  from  conventional  designs  employing  stacked  metal 
guard  rings  with  precision  insulating  spacers.  The  fabrication  costs  should,  in  general, 
be  less.  The  good  mobility  values  and  theoretical  fits  of  the  spectra  suggest  that  the 
exposed  insulators  between  the  coated  guard  rings  do  not  significantly  affect  the 
determination  of  ion  mobility  for  atmospheric  pressure  conditions. 


1.5  DC  Electric  Field  Effects  During  Measurement  of  Mono- 
polar  Charge  Density  and  Net  Space  Charge  Density 

This  NIST  study  briefly  examined  errors  that  arise  due  to  the  influence  of  a dc 
power  line  electric  field  during  measurements  of  monopolar  charge  densities  with 
an  aspirator-type  ion  counter  [9]  and  during  the  measurement  of  net  space  charge 
density  using  Faraday  cages  and  filters  [10,  11].  Where  appropriate,  measurements 
were  considered  with  the  instrumentation  in  and  above  the  ground  plane.  The  study 
made  use  of  theory  developed  by  Swann  for  terrestrial  electric  fields  and  air  ions  [12], 
and  recent  experimental  results  which  are  supportive  of  Swann’s  theory  even  for  cases 
where  there  is  significant  distortion  of  the  electric  field  due  to  space  charge  [13, 14];  the 
latter  case  can  be  found  near  high  voltage  dc  transmission  lines  [15,  16].  The  focus  of 
the  study  was  on  the  influence  of  the  electric  field.  Errors  which  can  occur  due  to  other 
factors  such  as  wind.  Coulomb  repulsion,  uncertainties  in  the  measurement  of  ion 
current  or  volumetric  air  flow  were  not  considered.  What  is  presented  in  the  following 
sections  are  excerpts  from  the  study  that  focuses  on  above-ground  measurements  of 
monopolar  charge  density  with  ion  counters  and  measurements  of  net  space  charge 
density  with  filters. 


1.6  Monopolar  Charge  Density  Measurements 
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Figure  6.  Schematic  view  of  aspirator-type  ion  coimters:  (a)  cylindrical  geometry, 
(b)  parallel-plate  geometry.  The  polarizing  potential  provided  by  the  battery  permits 
measurement  of  positive  charge  densities. 


1.6  Monopolar  Charge  Density  Measurements 


Schematic  views  of  cylindrical  and  parallel  plate  aspiratory-type  ion  counters  are 
shown  in  figure  6.  Under  ideal  conditions,  the  measured  monopolar  charge  density. 
Pm,  is  given  by  the  expression  [9] 


Pm  = I/M  (4) 

where  / is  the  measured  ion  current  to  the  collecting  electrode(s)  and  M is  the  laminar 
volumetric  flow  rate  of  the  air  through  the  ion  counter  due  to  the  action  of  the  fan. 
Equation  (4)  is  valid  for  ions  with  mobilities  equal  to  or  greater  than  the  “critical” 
mobility,  kc,  given  by  [9] 

k,  = Mco/CV  (5) 

where  C is  the  capacitance  of  the  electrodes,  Co  is  the  permittivity  of  free  space,  and 
V is  the  voltage  across  the  electrodes. 

It  can  be  shown  from  ion  trajectory  calculations  that  if  there  exists  an  electric  field 
external  to  the  ion  counter  that  transports  ions  of  a given  polarity  toward  the  ion 
counter,  the  number  of  ions  arriving  at  the  entrance  of  the  ion  counter  can  exceed  that 
due  to  the  aspirating  action  alone  by  a factor  of  two  or  more  [2] , It  can  also  be  shown 
that  if  there  is  adequate  electric  field  shielding,  i.e.,  the  collecting  electrode(s)  in  the 
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ion  counter  is  shielded  against  the  external  field,  the  excess  ions  will  travel  to  the 
housing  of  the  ion  counter  and  not  be  measured.  With  this  background  information, 
the  measurement  of  monopolar  charge  density  is  considered  for  an  ion  counter  above 
the  groimd  plane  in  the  presence  of  an  external  electric  field.  It  is  shown  that  if  the  ion 
counter  housing  is  appropriately  oriented  and  electrically  biased,  the  measurements 
can  be  made  free  of  errors  due  to  the  presence  of  the  electric  field. 

Figure  7 schematically  shows  ion  counters  in  two  orientations  above  the  ground  plane 
for  measuring  negative  charge  densities  in  the  presence  of  a uniform  vertical  electric 
field,  Eo-  The  housing  of  the  ion  counter  is  shown  electrically  floating  in  figures  7(a) 
and  7(d),  at  ground  potential  in  figures  7(b)  and  7(e)  and  biased  to  attract  negative 
ions  in  figures  7(c)  and  7(f).  The  number  of  induced  negative  charges  on  the  housing 
in  figures  7(b)  and  7(e)  depends  on  the  difference  in  potential  between  the  ambient 
space  potential  and  the  ion  counter  housing;  the  larger  this  difference,  which  is  a 
function  of  height  above  the  ground  plane,  the  fewer  are  the  number  of  negative 
charges. 

Electrically  floating  the  ion  counter  housing  at  ambient  space  potential  minimizes  the 
perturbation  of  the  electric  field.  However,  the  induced  negative  charges  near  the  ion 
counter  aperture  for  the  horizontal  orientation  [fig.  7(a)]  will  repel  some  approaching 
negative  ions  and  possibly  prevent  their  entry  into  the  ion  counter  [17].  When  the 
ion  counter  housing  is  at  ground  potential  or  has  an  attractive  potential  apphed  to  it 
[figs.  7(b)  and  7(c)],  the  number  of  approaching  negative  ions  repelled  by  the  induced 
negative  charges  is  reduced  or  eliminated. 

Figure  8 shows  experimental  results  demonstrating  ion  losses  by  the  mechanism  dis- 
cussed above  when  the  ion  counter  is  oriented  horizontally  about  0.4  m above  the 
ground  plane  and  beneath  a thin  horizontal  wire  2.1  m high  that  is  in  negative  corona. 
Measurements  of  negative  charge  densities  at  two  flow  rates  are  shown  as  a function 
of  voltage  apphed  to  the  ion-counter  housing  [13,  18].  During  the  measurements,  the 
unperturbed  electric  field  strength  in  the  ground  plane  directly  beneath  the  wire  was 
about  27  kV /m.  Further  details  of  the  experimental  arrangement  can  be  found  in  the 
references  that  are  cited. 

When  the  potential  of  the  housing  is  zero  or  positive  and  attracts  negative  ions,  the 
measured  charge  densities  are  nominally  constant  at  both  air  flow  rates.  The  scatter 
of  the  data  is  due,  in  part,  to  fluctuating  air  currents  in  the  laboratory  which  can  affect 
the  flow  of  ions  from  the  monopolar  line  to  the  ion  counter.  As  the  potential  of  the 
ion  counter  housing  approaches  the  ambient  unperturbed  space  potential  (determined 
with  an  electrostatic  voltmeter  and  radioactively  equilibrated  probe),  some  negative 
ions  are  repelled  by  the  negative  charges  on  the  housing  and  are  not  measured  as 
discussed  above.  It  is  noted  that  similar  experimental  results  have  been  obtained 
during  conductivity  measurements  aboard  airplanes  [19,  20]. 


1.6  Monopolar  Charge  Density  Measurements 
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Figure  7.  Schematic  views  of  ion  counter  above  ground  plane  in  vertical  and  horizontal 
orientations.  Induced  charges  produced  by  a uniform  vertical  electric  field,  Eq,  are 
shown  when  the  ion  counter  housing  is  at  space  potential  in  (a)  and  (d),  referenced 
to  ground  potential  in  (b)  and  (e)  [the  number  of  negative  charges  depends  on  the 
difference  in  potential  between  the  ion  counter  housing  and  ambient  space  potential] , 
and  biased  to  attract  negative  ions  in  (c)  and  (f). 
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Figure  8.  Above-ground  ion  density  measurements  at  two  air  flow  rates  and  as  a 
function  of  ion  counter  potential.  The  ion  counter  is  oriented  horizontally  and  the 
unperturbed  space  potential  at  the  height  of  the  aperture  is  near  — llkV. 

When  the  ion  counter  housing  is  oriented  vertically  and  electrically  biased  as  shown 
in  figures  7(d),  7(e),  and  7(f),  the  loss  of  approaching  negative  ions  due  to  repulsion 
from  like  charges  on  the  ion-counter  housing  is  expected  to  be  small  or  eliminated. 

It  appears  that  the  optimum  arrangement  for  above-ground  charge  density  measure- 
ments is  to  orient  the  ion  counter  horizontally  to  reduce  the  effects  of  rain  or  snow, 
and  to  apply  an  attractive  potential  to  the  ion  counter  housing  to  eliminate  ion  losses 
due  to  repulsion  by  like  charges  on  the  housing,  e.g.,  figure  7(c).  The  grounded  verti- 
cal geometry  shown  in  figure  7(e)  also  would  be  a useful  arrangement  in  the  absence 
of  adverse  weather  conditions.  It  is  noted  once  more  that  although  the  electric  field 
can  increase  the  number  of  ions  reaching  the  ion  counter  aperture  to  a value  in  ex- 
cess of  that  due  to  the  aspirating  action  of  the  ion  counter  alone,  the  extra  ions  will 
not  be  measured  provided  there  is  adequate  electric  field  shielding  of  the  collecting 
electrodes. 


1.7  Net  Space-Charge  Density  Measurements 

One  measurement  technique  for  determining  net  space  charge  density  is  known  as 
the  filtration  method  [10]  and  is  illustrated  in  figure  9.  Here,  an  air  sample  is  drawn 


1.7  Net  Space-Charge  Density  Measurements 
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Figure  9.  Schematic  view  of  filter  for  measuring  net  space  charge  density.  Both 
prefilter  and  HEP  A filters  are  isolated  from  ground  and  the  filter  housing.  To  prolong 
the  hfe  of  the  HEPA  filter,  it  is  common  practice  to  employ  a prefilter  made  of  such 
materials  as  stainless  steel  wool  and  glass  fibers. 

into  a metal  enclosure  with  a filter  through  which  all  the  air  passes.  The  filter  is 
electrically  isolated  from  the  enclosure  and  ground  so  that  the  current,  I,  resulting 
from  the  charge  removed  from  the  air  stream  by  the  filter  can  be  measured.  From 
this  current,  the  net  space  charge  density  can  be  calculated  from  the  expression 

Ps  = I/Mo  (6) 

where  M©  is  the  volumetric  air  flow  through  the  Alter.  This  method  of  measuring 
space  charge  density  assumes  that  the  filter  collects  almost  all  of  the  charge  in  the 
air  stream.  High-efficiency  air-particulate  (HEPA)  filters  are  adequate  for  this  apph- 
cation  because  tests  with  such  filters  indicate  that  the  collection  efficiency  for  small 
air  ions  can  exceed  99.8%  [21,  22]. 

Equation  6,  which  is  used  to  determine  net  space  charge  density,  is  derived  assuming 
that  ions  of  both  polarity  are  transported  into  the  filter  by  the  air  drawn  through  the 
filter.  However,  with  the  presence  of  an  external  electric  field,  ions  are  transported 
to  the  filter  aperture  under  the  combined  effects  of  the  air  motion  and  electric  field. 
It  is  noteworthy  that  the  electric  field  will  discriminate  against  the  arrival  of  ions  of 
one  polarity  and  favor  ions  of  the  opposite  polarity. 

Use  can  be  made  of  figure  7 for  identifying  an  optimum  arrangement  for  filter  measure- 
ments above  the  groimd  plane.  Taking  the  ion  counter  housing  in  figure  7 to  be  that 
of  a filter,  it  is  apparent  that  the  optimum  configuration  for  net  space  charge  density 
measurements  (but  not  monopolar  charge  density  measurements)  is  the  horizontal 
orientation  with  the  housing  electrically  floating  at  space  potential,  i.e.,  figure  7(a). 
For  this  configuration  ions  of  both  polarities  will  be  discriminated  against  by  the 
induced  charges  near  the  aperture.  If  the  positive  and  negative  ions  have  compara- 
ble mobihties,  the  percentage  of  positive  and  negative  ions  repulsed  by  the  induced 
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charges  will  also  be  comparable.  However,  even  with  equal  discrimination  against 
ions  of  both  polarities,  the  net  space  charge  density  measurement  can  be  adversely 
affected  unless  the  net  space  charge  is  zero. 

Measurements  of  net  space  charge  density  using  a filter  above  the  ground  plane,  with 
the  filter  housing  at  space  potential,  is  reported  by  Crozier  [23].  However,  details 
of  the  experimental  arrangement  are  not  provided.  Extra  ions  which  arrive  at  the 
aperture  because  of  the  electric  field  will  not  be  measured  if  the  filter  element  is 
adequately  shielded  by  the  filter  housing.  The  electrical  shielding  can  be  checked  by 
noting  whether  there  is  any  ion  current  to  the  filter  when  Mo  is  zero. 


1.8  Conclusion 

This  study,  described  in  more  detail  elsewhere  [2],  has  shown  that  the  influence  of 
a dc  electric  field  on  measurement  of  monopolar  space  charge  densities  with  an  ion 
counter  can  be  made  neghgible  when  the  measurements  are  performed  in  the  ground 
plane,  and  oriented  vertically  or  horizontally  above  the  ground  plane  with  an  attrac- 
tive potential  apphed  to  the  ion  counter  housing.  This  assumes  that  the  collecting 
electrode(s)  in  the  ion  counter  is  adequately  shielded  from  the  external  electric  field. 

Measurements  of  net  space  charge  density  with  a Faraday  cage  (not  discussed  in 
this  report)  appear  to  be  significantly  influenced  by  a large  field  except  when  the 
space  charge  is  predominantly  monopolar  or  when  composed  of  large  ions  (aerosols) 
with  small  mobilities.  The  optimum  configuration  for  the  filter  is  an  above-ground 
horizontal  orientation  with  the  filter  housing  at  ambient  space  potential. 
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2.1  Introduction 

The  objectives  of  this  project  are  the  development  of  measurement  capabilities  and 
providing  fundamental  data  as  part  of  the  Department  of  Energy’s  basic  research  con- 
cerned with  the  development  and  evaluation  of  advanced  compressed-gas-insulation 
technology. 

To  reduce  space  requirements  and  improve  the  efficiency  of  high-voltage  transmis- 
sion systems,  the  electric  power  industry  has  turned  toward  more  extensive  use  of 
compressed-gas  insulation.  To  design  meaningful  tests  of  system  performance  and 
establish  specifications  for  the  quahty  of  materials  used  in  such  systems,  more  in- 
formation is  needed  about  the  fundamental  physical  and  chemical  processes  which 
lead  to  insulation  deterioration  and  electrical  breakdown.  The  research  includes  ap- 
plications of  gas  chromatography-mass  spectrometry  to  characterize  corona  discharge 
by-products;  and  the  acquisition  of  fundamental  data,  e.g.,  reaction-rate  coefficients, 
corona-inception  voltages,  production  rates  of  corona  by-products,  the  effects  of  con- 
taminants on  discharge  initiation,  and  the  rates  of  discharge-induced  decomposition 
of  the  gas. 

This  report  highlights  two  significant  technical  activities.  The  first  is  concerned  with 
measurement  of  colhsional  electron-detachment  and  ion-conversion  cross  sections  for 
the  interaction  of  the  negative  ions  SFg  , SFj  , and  F“  with  SFe  molecules.  The  mea- 
sured cross  sections  were  used  to  calculate  electron-detachment  and  ion-conversion 
reaction  rate  coefficients  which  were  in  turn  used  in  a model  to  interpret  previous 
results  of  negative  ion  measurements  in  imiform-field  drift  tubes  and  investigations 
of  electrical  discharge  initiation  probabihties  at  relatively  high  pressures  in  SFe  insu- 
lated power  systems.  The  experimental  part  of  this  work  was  conducted  in  a collab- 
oration with  scientists  in  the  Department  of  Physics,  College  of  William  and  Mary, 
Williamsburg,  Virginia.  The  results  of  this  work  have  been  accepted  for  publication 
in  two  archival  papers  [24,  25]  and  have  also  been  presented  in  three  recent  conference 
publications  [26-28]. 

The  second  major  technical  activity  to  be  highlighted  concerns  the  development  of 
a new  technique  to  measure  the  stochastic  properties  of  pulsating  corona  or  partial 
discharges.  This  technique  is  unique  in  that  it  allows  direct  measurement  of  various 
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conditional  pulse-height  and  pulse-time  interval  distributions  that  enable  a quanti- 
tative determination  of  the  degree  of  correlation  among  successive  partial  discharge 
pulses.  A detailed  description  of  the  measurement  system  has  been  accepted  for 
archival  publication  [29].  Results  from  apphcation  of  this  method  to  an  investigation 
of  the  stochastic  behavior  of  negative  corona  (Trichel)  pulses  in  electronegative  gas 
mixtures  have  been  presented  at  three  different  technical  conferences  [30-32]. 


2.2  Collisional  Electron-Detachment  and  Ion- Conversion 
Processes  in  SFg 

2. 2. 1 Introduction 

It  has  long  been  hypothesized  that  colhsional  electron  detachment  from  negative  ions 
may  be  the  cause  of  discharge  initiation  in  SFe-insulated  high-voltage  equipment. 
However,  uncertainties  continue  to  exist  concerning  the  identities  of  the  significant 
negative  ions,  and  the  magnitudes  of  the  detachment  cross  sections  and  correspond- 
ing reaction  coefficients.  Discharge-inception  studies,  which  usually  assume  that  SFg 
is  the  negative  ion  of  interest,  have  predicted  colhsional-detachment  cross  sections 
of  SF0  in  SFe  with  thresholds  ranging  from  1 to  8eV  [33].  Colhsional-detachment 
coefficients  determined  from  low-pressure  drift-tube  experiments  [34,  35]  lie  several 
orders  of  magnitude  above  those  calculated  from  breakdown  data  [33],  while  ex- 
hibiting an  unexplained  pressure  dependence.  Additionally,  evidence  exists  that  the 
destruction  of  negative  ions  in  SFe  is  dominated  by  colhsional  dissociation  into  ionic 
fragments  [34,  36-40]  or  by  charge- transfer  processes  rather  than  by  the  loss  of  an 
electron.  However,  the  exact  nature  of  these  ion-conversion  processes  have  not  pre- 
viously been  determined. 

In  this  report  we  summarize  results  from  the  first  direct  measurements  of  absolute 
cross  sections  for  electron-detachment  and  ion-conversion  processes  involving  inter- 
actions of  SFg  , SF5  , and  F~,  with  SFe.  These  cross  sections  are  used  to  calculate 
electron-detachment  and  ion-conversion  reaction  coefficients  as  functions  of  electric 
field-to-gas  density  ratios  (E/N)  for  the  reactions  listed  in  Table  1.  We  then  discuss 
the  relevance  of  these  results  to  the  interpretation  of  data  from  uniform-field  drift- 
tube  measurements  and  measurements  of  electrical-discharge  initiation  processes. 


2.2.2  Cross  Sections 

Details  of  the  experimental  method  used  to  measure  the  cross  sections  reported  here 
are  given  elsewhere  [25,  41],  and  will  not  be  covered  here.  We  shah  focus  here  only 
on  the  results  obtained  for  colhsions  of  SFg  , SF5  , and  F”  on  the  SFe  target  gas. 
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Table  1.  Collisional  processes  for  whicii  cross  sections  have  been  presented  in  the 
present  work 


cross  section  (crj) 

reaction 

SFg  -f-  SFe  — > e”  -f  SFe  + SFe 

SFg  -f  SFe  — > e~  -j-  SF5  -f  SFe 

0-3 

F~  -1-  SFe  ->  e“  + F -1-  SFe 

04 

SFg  + SFe  F-  + SF5  + SFe 

^5 

SFg  -|-  SFe  — > F~  + SF4  -f  SFe 

SFg  + SFe  ->  SFg-  -b  F -f  SFe 

07 

SFg  + SFe  — > charge-transfer  products  -f  SFe 

(^8 

SFg  -f-  SFe  — > charge-transfer  products  + SF5 

(Tg 

F”  + SFe  — ♦ F + SFe  (see  ref.  35) 

although  measurements  were  also  made  using  other  target  gases  [25],  namely  He,  Ne, 
Ar,  Kr,  and  Xe. 

The  experimentally  determined  center-of-mass  energy  dependencies  of  the  colhsional 
electron-detachment  cross  sections,  <T,(€cin)  for  SFg  , SF5  , and  F"  on  an  SFe  target  are 
shown  in  Fig.  10.  The  thresholds  for  the  “prompt”  colhsional  electron-detachment 
processes 

SFg  -f  SFe  — > neutral  products  -t-  e“ 

SF5  + SFe  — > neutral  products  + e” 

as  indicated  by  the  vertical  arrows  in  Fig.  10,  are  seen  to  be  quite  high  (approxi- 
mately 90 eV  in  both  cases),  while  F~  has  a lower  detachment  threshold  of  about 
8eV.  At  center-of-mass  colhsion  energies  below  the  indicated  thresholds  there  is  no 
evidence  to  suggest  finite  (non  zero)  values  for  the  <Ti(ecin)’s.  Similar  high  values  for 
the  colhsional  electron-detachment  thresholds  were  found  for  colhsions  of  SFg  and 
SFg  ions  with  rare  gas  targets  [25].  The  implication  of  these  results  as  shown  below 
is  that  prompt  electron  detachment  from  either  SFg  or  SFg  cannot  be  important  for 
electrical  discharge  conditions  because,  even  for  the  highest  E/N  values  which  could 
conceivably  occur  in  a gas  discharge  gap,  only  an  insignificant  fraction  of  the  negative 
ions  could  acquire  a kinetic  energy  of  90  eV  or  more. 

Cross  sections  for  the  ion-conversion  processes  (reactions  (4)-(9))  listed  in  Table  1 
are  shown  in  Fig.  11  together  with  extrapolations  down  to  the  thresholds  which 
are  used  in  the  calculation  of  corresponding  rate  coefficients.  The  bases  for  these 
extrapolations  are  discussed  later.  The  solid  symbols  indicate  cross  sections  due 
to  colhsional-induced-dissociation  (CID)  processes,  and  the  open  symbols  indicate 
cross  sections  for  charge- transfer  processes  (including  charge-transfer  decomposition). 
There  is  evidence  that  the  charge-transfer  process  involving  SFg  on  SFe  is  predomi- 
nately dissociative  [42,  43],  and  that  charge  transfer  involving  F“  on  SFe  may  lead  to 


ELECTRON  DETACHMENT 


20 


2.  GASEO US  DIELECTRICS  RESEARCH 


Figure  10.  Collisional  electron-detachment  cross  sections  for  F , SF5  , and  SFg  on 
SFe  target  gas  as  a function  of  center-of-mass  energy. 
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Figure  11.  Measured  cross  sections  for  collision-induced  ion-conversion  processes  in 
SFe  target  gas:  (o)  F“  from  SFj ; (□)  F"  from  SFg  ; (A)SF5  from  SFg  ; (o)  ions  due 
to  charge-transfer  reactions  of  SF5  ; (□)  ions  due  to  charge-transfer  reactions  of  SFg  ; 
and  (v)ioiis  due  to  charge- transfer  reactions  of  F".  The  two  exponentially-decaying 
curves  on  the  left  side  of  the  figure  represent  the  kinetic-energy  distribution  of  Eq.  (11) 
for  7 = 1.0  and  7 = 0.5  scaled  relative  to  each  other. 
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both  SFg  and  SF5  [34,  42].  It  is  interesting  to  note  that  the  ion-conversion  processes 
predominate  for  Ccm  < 200eV,  i.e.,  upon  colhsion  of  an  SFg  or  SFg  ion  with  SFe,  ion 
conversion  is  much  more  likely  than  electron  detachment. 

2.2.S  Reaction  Rates  and  Coefficients 

The  analysis  of  rates  for  chemical  processes  in  drift  tubes  or  electrical  discharges 
requires  expressing  inelastic-collision  probabihties  in  terms  of  rate  coefficients  rather 
than  cross  sections.  For  a process  where  the  projectiles  have  a velocity  distribution, 
/(u),  the  rate  coefficient,  k,  is  given  by  an  averaging  of  the  colhsion  cross  section  over 
all  velocities,  namely 

roo 

k=  a{v)vf{v)dv, 

Jo 

where  a(u)  is  the  velocity-dependent  cross  section  for  the  process  and  f{v)  is  subject 
to  the  normalization 

/ f{v)dv  = 1.  (8) 

Jo 

For  measurements  where  charged  particles  are  accelerated  in  a gas  by  an  electric 
field  (such  as  in  a drift  tube),  the  reaction  coefficient,  «,  is  defined  as  the  reaction 
probabihty  per  unit  length  in  the  direction  of  the  electric  field,  and  is  related  to  the 
rate  coefficient  by 

K k 
N ^ 

where  N is  the  target-gas  number  density  and  Vd  is  the  charged  particle  drift  velocity. 
For  the  special  case  of  colhsional  electron  detachment  from  a negative  ion,  the  reaction 
coefficient  is  referred  to  as  the  detachment  coefficient,  6 (i.e.  ac/N  = 6/N  for  colhsional 
electron  detachment). 

Reaction  coefficients  for  detachment  and  ion-conversion  processes  involving  SF0  , SF5  , 
and  F“  in  SFe  can  be  derived  from  the  measured  cross  sections,  cri(ecni)5  presented  in 
Figs.  10  and  11  using 

AT*  1 

^ = / o-.(€L)/(cL)<^eL,  (10) 

N mvd  Jo 

where  cl  is  the  projectile  energy  in  the  lab  frame,  m is  the  mass  of  the  negative 
ion,  and  /(cl)  and  cr{ei,)  are  the  ion  kinetic-energy  distribution  and  the  process  cross 
sections,  respectively.  The  ion  drift  velocity  is  assiuned  to  be  given  by  Vd  = /aE,  where 
fi  is  the  ion  mobility  and  E is  the  apphed  electric-field  strength. 

It  should  be  noted  that  the  determination  of  reaction  coefficients  from  colhsional  cross 
sections  suffers  from  certain  difficulties,  the  greatest  arising  from  the  assumed  form 
of  the  ion  kinetic-energy  distribution,  /(cl)  [44].  While  determination  of  ion  kinetic- 
energy  distributions  has  received  considerable  theoretical  [45]  and  experimental  [46] 
attention,  the  theoretical  work  is  hampered  by  the  lack  of  detailed  ion-atom  (or 
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molecule)  potential-energy  surfaces  and  the  experimental  work  suffers  from  a lack  of 
reliability  [44].  Accurate  direct  measurements  of  ion-velocity  distributions  have  been 
demonstrated  in  recent  optical-probing  experiments,  [47]  but  to  date  no  experimental 
data  are  available  for  the  kinetic-energy  distributions  of  SFe  , SFj  , or  F“  in  SFe. 


In  general,  experimental  work  has  indicated  that  ions  with  masses  less  than  or  equal 
to  that  of  the  molecules  of  the  gas  in  which  they  are  moving  exhibit  kinetic-energy 
distributions  with  high-energy  tails  [44,  48].  Differences  in  the  velocity  distributions 
at  high  energies  in  Eq.  (10)  will  obviously  have  a large  effect  upon  the  calculated 
reaction  coefficients  derived  from  cross  sections  with  threshold  energies  considerably 
in  excess  of  the  average  ion  kinetic  energies. 

An  example  of  the  large  differences  in  calculated  values  of  reaction  coefficients  which 
can  occur  when  different  energy  distributions  are  assiuned  is  shown  in  Fig.  12,  where 
the  colhsional-detachment  coefficient  for  F“  in  SFe  ((^s)  has  been  calculated  as  a func- 
tion of  E/N  using  the  cross  section  data  (era)  from  Fig.  10  and  several  different 
indicated  energy  distributions.  The  solid  lines  in  Fig.  12  represents  detachment  coef- 
ficients calculated  using  the  kinetic-energy  distribution  of  Kagan  and  Perel  [33,  49], 


/(^l)  = 


TTVd 


exp 


■KTriv^J  ’ 


(11) 


which  assmnes  that  charge  exchange  is  the  dominant  ion-molecule  interaction.  For 
the  standard  Kagan  and  Perel  distribution,  7 = 1.0.  However,  as  will  be  discussed 
later,  better  agreement  between  ion-conversion  reaction  coefficients  calculated  here 
and  those  from  analysis  of  drift-tube  results  is  obtained  by  assuming  7 = 0.5  which 
introduces  a larger  high-energy  tail  in  the  distribution. 


The  dashed  line  in  Fig.  12  was  obtained  using  a Maxwellian  speed  distribution  [33] 
of  the  form. 


/(«l)  = 3 


/3m  €l  / 

'-3eL\ 

/ . ( 

< 2€  J 

(12) 


where  the  mean  energy  of  the  ion  in  the  lab  frame  is  given  by  [50] 


ZkT  m 2 Mo 

and  M is  the  mass  of  the  coUision-gas  molecules.  This  distribution  has  been  used 
previously  when  analyzing  discharge-inception  data  [33]  and  is  similar  to  a strongly 
anisotropic  velocity  distribution  derived  by  Skullerud  [51]  for  ions  drifting  in  a gas 
composed  of  molecules  of  the  same  mass  as  the  ions  under  high  electric  fields. 

In  addition  to  the  uncertainty  associated  with  the  assumed  distribution  function,  an- 
other source  of  uncertainty  in  deriving  reaction  coefficients  from  cross-section  data 
is  the  choice  of  experimentally  determined  values  of  ion  mobilities  (/x).  Several  pre- 
vious mobility  measurements  [52-55]  for  SFg  and  SF5  in  SFe  are  not  in  complete 
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E/N  (10-21  Vm2) 


Figure  12.  Calculated  collisional  electron-detachment  coefficients  for  F“  on  SFe  gas 
using  Eq.  (10)  and  and  assuming  different  ion  kinetic-energy  distributions: 

( ) Kagan  and  Perel;  ( ) Maxwellian. 
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agreement.  However,  uncertainties  in  calculated  reaction  coefficients  using  different 
mobilities  are  significantly  smaller  than  the  uncertainties  due  to  the  use  of  different 
energy  distributions  as  discussed  above.  For  the  remainder  of  this  paper,  the  mobil- 
ities used  for  SFg  and  SFJ  are  those  reported  by  Brand  and  Jungblut,  [52]  and  for 
F",  the  values  of  Nakamura  [55]. 

As  expected,  the  extremely  high  apparent  thresholds  for  prompt  colhsional  electron 
detachment  from  SF^  and  SF5  yield  detachment  coefficients  from  Eq.  (10)  that  are 
insignificantly  small  in  the  E/N  range  of  interest  here.  Simple  estimates  of  these  co- 
efficients using  a Kagan  and  Perel  velocity  distribution  in  Eq.  (10)  indicate  that  the 
detachment  coefficients  for  SFg  and  SF5  will  be  tens  of  orders  of  magnitude  below 
the  detachment  coefficients  determined  in  drift-tube  experiments  and  the  detach- 
ment coefficients  calculated  for  F“  (Fig.  12).  Thus  one  must  conclude  that  prompt 
colhsional  electron-detachment  processes  for  SFg  (and  SFj ) in  SFe  cannot  be  signifi- 
cant reactions  for  production  of  electrons  in  discharge-inception  processes,  if,  in  fact, 
cri(ecm)  is  zero  below  the  apparent  thresholds.  It  should  be  noted  that  this  conclusion 
is  independent  of  the  assiuned  kinetic-energy  distribution  or  the  ion-mobility  data 
used. 

It  is  clear  from  Fig.  12,  however,  that  colhsional  electron  detachment  from  F"  in  SFe 
is  a significant  process,  thus  suggesting  that  previously  observed  electron-detachment 
processes  due  to  motion  of  negative  ions  in  SFe  are  very  likely  from  F“.  This  agrees 
with  earher  work  [56]  which  indicated  that  detachment  in  SFe  was  predominately  from 
F“,  and  also  with  recent  reanalysis  of  pulsed-electron  avalanche  experiments  [57].  In 
fact,  the  observed  threshold  for  electron  detachment  from  F“  in  SFe  near  8eV  is 
consistent  with  the  hypothesized  thresholds  predicted  by  discharge-inception  experi- 
ments [33]. 

Other  discharge-inception  experiments  have,  however,  suggested  that  negative  clus- 
ter ions  [58]  involving  H2O  or  HF  may  be  responsible  for  discharge  initiation  in 
SFe.  In  fact,  recent  measurements  [59]  of  negative  ions  produced  in  negative  corona 
discharges,  have  indicated  the  presence  of  several  types  of  cluster  ions  of  the  forms 
F~(HF)„,  0H“(H20)n,  and  SFg  (HF).  Further  investigations  of  the  colhsional  detach- 
ment cross  sections  for  these  ions  are  necessary  to  determine  their  role  in  discharge 
initiation. 

The  conclusions  drawn  above  depend  upon  the  assumption  that  = 0 at  energies 
below  the  apparent  threshold  in  Fig.  10.  If  one  assimies  that  cr{€i)  = 0.1  (i.e., 

the  experimental  uncertainty)  for  energies  which  extend  down  to  the  thermodynamic 
threshold  for  electron  detachment,  then  detachment  coefficients  for  SF^  and  SFJ 
are  found  to  be  of  the  same  order  of  magnitude  as  those  determined  by  drift-tube 
experiments.  However,  detachment  coefficients  derived  with  such  an  assumption  are 
not  compatible  with  previously  observed  pressure  dependencies  observed  in  drift-tube 
experiments  [34,  35]  as  discussed  in  the  next  section. 
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In  order  to  calculate  the  reaction  coefficients  for  the  ion-conversion  processes  listed  in 
Table  1,  it  is  necessary  to  extrapolate  the  measured  cross  sections  down  to  assiuned 
thresholds  at  lower  energies.  The  extrapolations  used  for  the  subsequent  calcula- 
tions are  shown  in  Fig.  11.  These  extrapolations  were  chosen  to  agree  with  known 
thermodynamic  thresholds  and  to  minimize  the  discrepancies  with  previously  deter- 
mined reaction  coefficients  as  discussed  below.  The  thresholds  for  production  of  F“ 
from  SFg  and  SFj  (<74  and  <75)  were  determined  to  be  2.0 eV  and  1.5 eV,  respectively, 
by  using  the  observed  thresholds  for  F“  production  from  colhsions  of  SFg  and  SF5 
with  the  rare-gas  targets  [25].  The  cross  sections  for  SFj  production  (as)  from  SFg 
were  extrapolated  down  to  the  thermodynamic  threshold  of  1.35  eV.  For  the  charge- 
transfer  reaction  involving  F~  and  SFe,  the  cross  sections  were  extrapolated  down 
to  the  thermodynamic  threshold  of  2.25  eV  under  the  assumption  that  the  primary 
product  is  SFg  . The  other  charge- transfer  cross  sections  (<77  and  ag)  were  both  ex- 
trapolated down  to  a threshold  near  3eV  which  corresponds  to  the  thermodynamic 
threshold  for  a symmetric  charge  transfer  between  SFg  and  SFe  as  suggested  by 
Hay  [60].  There  is  a large  uncertainty  in  these  last  assumed  thresholds  since  the 
identity  of  the  charge-transfer  products  in  these  processes  are  indistinguishable  in 
the  present  experiment. 

The  calculated  reaction  coefficients  for  processes  5,  6 and  9 of  Table  1 are  shown  in 
Figs.  13,  14,  and  15  respectively,  along  with  reaction  coefficients  for  the  same  reac- 
tions as  determined  by  previous  drift-tube  experiments.  The  solid  lines  represent  the 
coefficients  calculated  using  the  standard  Kagan  and  Perel  distribution  (i.e.  7 = 1) 
shown  in  Eq.  (11).  Note  that  these  calculated  reaction  coefficients  all  fall  substantially 
below  those  determined  previously  despite  the  fact  that  the  extrapolated  thresholds 
for  these  cross  sections  were  all  assumed  to  be  thermodynamic  thresholds.  Any  rea- 
sonable change  in  the  assumptions  concerning  the  reaction  thresholds  or  the  behavior 
of  the  cross  sections  near  threshold  would  necessarily  cause  the  reaction  coefficients  to 
be  even  smaller,  thus  implying  that  the  discrepancies  cannot  be  resolved  by  changing 
the  assumed  cross-section  thresholds  or  extrapolations. 

Despite  the  fact  that  the  Kagan  and  Perel  distribution  produces  the  largest  coeffi- 
cients of  any  of  the  commonly  used  ion  kinetic-energy  distributions,  better  agreement 
can  be  obtained  between  our  calculated  coefficients  and  those  from  previous  exper- 
iments if  one  assumes  that  the  kinetic-energy  distribution  has  a longer  high-energy 
tail,  in  agreement  with  the  previous  discussion  of  kinetic-energy  distributions.  The 
Kagan  and  Perel  distribution  in  Eq.  (11)  can  be  conveniently  altered  by  allowing 
7 to  vary  between  0 and  1.  The  dashed  lines  in  Figs.  13  to  15  represent  reaction 
coefficients  calculated  with  7 = 0.5,  while  in  Fig.  11  the  relative  magnitudes  of  the 
two  distributions  (7  = 1 and  7 = 0.5)  are  shown  for  comparison. 

Obviously,  the  curves  in  Figs.  13-15  with  7 = 0.5  are  in  better  agreement  with  the 
previously  reported  coefficients  than  are  the  curves  calculated  using  7 = 1.  This  may 
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E/N  (10-21  Vm2) 


Figure  13.  Calculated  reaction  coefficients  for  the  reaction  SFg  -f-  SFe  — > SF5  -f 
F + SFe  using  measured  cross-section  data  and  a Kagan  and  Perel  energy 

distribution.  The  symbols  are  previously  calculated  reaction  coefficients  for  the  same 
process  derived  from  uniform-field  drift-tube  data. 
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Figure  14.  Calculated  reaction  coefficients  for  the  reaction  SF5  + SFg  — F“  + 
SF4  + SFe  using  measured  cross-section  data  <75(6l)  and  a Kagan  and  Perel  energy 
distribution.  The  symbols  are  previously  calculated  reaction  coefficients  for  the  same 
process  derived  from  uniform-field  drift-tube  data. 
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E/N  (10*21  V m2) 


Figure  15.  Calculated  reaction  coefficients  for  the  reaction  F~  + SFe  — > F + SFg 
using  measured  cross-section  data  cr9(cL)  and  a Kagan  and  Perel  ion  kinetic-energy 
distribution.  The  symbols  are  previously  calculated  reaction  coefficients  for  the  same 
process  derived  from  uniform-field  drift-tube  data. 
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indicate  that  previously  assumed  energy  distributions  need  to  be  modified  [33,  35]. 
However,  one  must  note  that  the  reaction  coefficients  derived  from  drift-tube  data 
are  model  dependent  and  that  only  reactions  5,  6,  and  9 (table  1)  are  assumed  in  the 
previous  analysis  of  data  from  drift  tubes.  Thus  discrepancies  may  also  arise  because 
this  commonly  used  model  does  not  consider  the  colhsion-induced  dissociation  of  SFg 
into  F“  -f  SFs  (reaction  4).  This  reaction  is  found  here  to  be  significant  (see  Fig.  11) 
and  its  omission  may  produce  errors  in  the  reaction  rates  derived  from  drift-tube 
data. 

The  calculated  reaction  coefficients  for  reaction  4 (and  for  reactions  7 and  8)  are 
shown  in  Fig.  16  using  Kagan  and  Perel  distributions  with  7 = 1.0  and  0.5.  As  stated 
before,  no  previously  determined  coefficients  for  these  processes  exist  for  comparison. 


2.2.4  Model  for  Electron  Detachment  from  Ion  Drift  in  SFq 


A different  interpretation  of  the  processes  which  lead  to  detachment  coefficients  de- 
rived from  drift-tube  data  [34,  35]  and  their  observed  pressure  dependence  can  be 
obtained  if  one  assumes  that  electron  production  in  a drift  tube  is  not  due  primarily 
to  prompt  electron  detachment  from  SFg  , but  arises  from  either  detachment  from  F“ 
or  from  colhsionally-excited,  energetically-unstable  (SFe  )*  ions  as  suggested  by  the 
previously  presented  data.  A model  can  then  be  developed  which  assiunes  that  under 
drift-tube  conditions,  prompt  detachment  from  SF5  and  SFg  is  insignificant,  and  that 
a steady  state  condition  exists  for  intermediate  products  (i.e.  ^[SFj ]/dt  = d[F”]/dt 
= d[{SFQ)*]/dt  = 0).  If  one  analyzes  the  drift-tube  data  as  done  previously  [34,  35], 
assiuning  that  all  ejected  electrons  come  from  SFg  ions,  then  the  measured  electron- 
production  rate  can  be  written  in  terms  of  an  “effective”  detachment  coefficient  (^eff) 
according  to  the  expression 

f = ..(^)n[SFJ1,  (14) 


where  [e]  is  the  number  of  free  electrons  produced  per  unit  volume  which  can  be 
detected,  Vd  refers  specifically  to  the  drift  velocity  of  SF^  in  SFe,  and  <5eff/N  is  de- 
termined by  analysis  of  drift-tube  data  [34,  35].  Assuming  that  the  simplified  set  of 
processes  indicated  in  Table  2 dominate  in  a drift  tube,  an  expression  for  ^eff/N,  can 
be  found  in  terms  of  the  relevant  rate  coefficients  by  solving  the  set  of  coupled  rate 
equations.  The  result  is 


Ka  = Vd 


+ k 


‘-10 


12 


ki2  + AtiiN 


(15) 


The  effective  detachment  coefficient  given  by  Eq.  (15)  is  seen  to  consist  of  two  terms, 
a pressure-independent  term  which  depends  upon  various  ion-conversion  and  direct- 
detachment  processes  involving  F“,  and  a pressure- dependent  term  which  depends 
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200  300  400  500  600 

E/N  (10*21  V m2) 


Figure  16.  Calculated  reaction  coefficients  using  a Kagan  and  Perel  ion  kinetic-energy 

distribution  for  the  following  reactions:  ( ) SFg  -f-  SFe  — > F“  -|-  SF5  -|-  SFe; 

( — • — ) SFg  -H  SFe  — > charge-transfer  products;  and  ( ) SF5  -I-  SFe  - — > 

charge- transfer  products. 
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Table  2.  Proposed  ion-molecule  reactions  for  drift-tubes  containing  SFr 


SFg  -b  SFe 

-^SF--f  F-f  SFe 

ke  ' 

SFg  -|-  SFe 

F-  + SFe  + SFe 

k^ 

> Dissociative  Ion  Conversion 

SFg  -f  SFe 

— >■  F -f  SF4  + SFe 

^5  . 

F-  + SFe - 

-F  + SFg- 

/jg 

Charge  Transfer 

F-  -f  SFe  - 

— > neutrals  -f  e” 

h 

e~  Detachment 

SFg  + SFe 

— .(SFg-)- -f  SFe 

^10 

Excitation 

(SFg-)-  -f  SFe  — > SFg  + SFe 

k\i 

De-excitation 

(SFe-)-  — 

SFe  + e 

ki2 

Auto-detachment 

upon  the  rates  for  collisional  relaxation,  excitation  and  auto-detachment  of  (SFg  )*. 
This  expression  is  more  complex  than  those  previously  derived  [33-35]  which  assume 
that  direct  detachment  from  SFg  was  the  sole  source  of  electrons. 

If  AriiN  is  approximately  the  colhsion  frequency  of  SFg  in  SFe  (e.g.,  ~ 10®/sec  at 
1 kPa)  and  is  on  the  order  of  the  inverse  of  the  excited-state  lifetime  (r  ~ lOfis  to 
2ms),  [61]  then  AthN  >>  ki2  and  the  model  predicts  an  inverse  pressure  dependence 
for  6efr/N  at  low  N.  This  inverse  pressure  dependence  is  consistent  with  previous  drift- 
tube  measurements  [34,  35].  O’Neill  and  Graggs  [34]  also  report  no  detachment  from 
F“  or  SFj  at  low  pressures  in  agreement  with  the  dominance  of  the  second  term  of 
Eq.  (15)  for  smaller  N.  The  model  proposed  here  is  consistent  with  the  observed  [35] 
variations  in  ^eff/N  with  the  “age”  of  the  SF^  ions  if  a substantial  fraxtion  of  SFe 
anions  are  initially  in  excited  energetically- unstable  autodetaching  states. 

At  higher  pressures,  the  first  term  on  the  right  side  of  Eq.  (15)  dominates,  giving 
a <5eff/N  that  is  essentially  pressure  independent.  In  this  pressure  regime,  electron 
production  involves  mainly  process  (3)  with  a threshold  of  8eV.  These  results  are 
consistent  with:  1)  the  lack  of  pressure  dependence  [33,  57]  for  6eff/N  suggested  from 
the  analysis  of  high-pressure  electrical  discharge  initiation  data,  and  2)  previously 
discussed  results  [56]  suggesting  detachment  in  SFe  is  predominately  from  F~. 

Ideally,  one  would  like  to  calculate  the  effective  detachment  coefficients  using  Eq.  (15) 
to  compare  with  the  previously  determined  drift-tube  measurements.  However,  values 
for  kioi  ku  and  k\2  are  not  available,  so  only  the  contribution  to  <5eff/N  from  the  first 
term  of  Eq.  (15)  can  be  calculated  using  the  reaction  coefficients  derived  above.  The 
solid  curve  in  Fig.  17  shows  the  contribution  from  the  first  term  of  Eq.  (15)  using 
a Kagan  and  Perel  distribution  with  7 = 1.0.  Note  that  the  magnitude  of  the  solid 
curve  is  similar  to  that  of  the  effective  detachment  coefficient  derived  from  discharge- 
inception  experiments  [33]  (dot-dashed  curve)  but  is  substantially  smaller  than  the 
coefficients  derived  from  drift-tube  experiments  [34,  35]  (symbols).  If  one  uses  the 
reaction  coefficients  derived  with  a Kagan  and  Perel  energy  distribution  for  7 = 0.5 
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E/N  (10-21  V m2) 


Figure  17.  Effective  detachment  coefficients  determined  from  drift-tube  experiments 
at  different  pressures  are  shown  as  identified  by  the  symbols  in  the  key.  Effective 
detachment  coefficients  calculated  from  discharge-inception  data  from  Ref.  [33]  are 
also  shown  ( — • — ).  The  contribution  to  the  effective  detachment  coefficients  from 
the  first  term  of  Eq.  (15)  are  shown  for  7 = 1.0  ( ) and  for  7 = 0.5  ( ). 
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(dashed  curve),  then  the  <5efr/N  derived  from  the  first  term  of  Eq.  (15)  becomes  of  the 
same  order  of  magnitude  as  the  coefficients  derived  from  the  highest  pressure  drift- 
tube  experiments.  The  fact  that  the  dashed  curve  actually  lies  above  the  smallest 
measured  drift-tube  values  indicates  that  taking  7 = 0.5  may  overestimate  the  high- 
energy  tail  for  the  kinetic-energy  distribution. 

It  is  also  important  to  note  that  the  above  model  gives  a reasonable  pressure  depen- 
dence for  ^eff/N  without  invoking  “stabilization”  of  SFg  via  (SF6)SF6  dimer  forma- 
tion as  proposed  by  Hansen  et  al.  [35]  Dimer  formation  is  irrelevant  here  because: 
1)  it  was  introduced  to  compete  with  the  direct  SFg  colHsional-detachment  process, 
which  is  of  doubtful  significance  and,  2)  as  can  be  shown,  it  cannot  yield  the  observed 
pressure  dependence.  There  are  also  data  which  suggest  that  (SF6)SF6  formation  is 
very  slow  [54]  and  decreases  with  increasing  E/N  [62]. 


2.2.5  Conclusions 

The  measured  cross  sections  for  colhsional  electron  detachment  and  ion  conversion  of 
negative  ions  in  SFg  have  been  used  in  a theoretical  model  which  invokes  detachment 
from  long-lived,  energetically-unstable  states  of  colHsionally  excited  SFg  to  explain 
the  pressure  dependence  of  previously  measured  detachment  coefficients  and  the  high 
apparent  detachment  thresholds  implied  by  analysis  of  breakdown-probability  data 
for  SFe.  The  model  suggests  that  measured  effective  detachment  coefficients  depend 
upon  many  different  reactions  rates,  thus  implying  that  detachment  processes  in  SFe 
are  more  complex  than  previously  assmned.  At  high  pressures,  measured  detachment 
coefficients  appear  to  depend  primarily  upon  the  rates  for  ion-conversion  and  direct- 
detachment  processes  involving  F“,  consistent  with  earher  suggestions. 


2.3  Stochastic  Properties  of  Partial  Discharges 

2.S.1  Introduction 

Pulsating  electrical  discharge  phenomena  such  as  corona  or  partial-discharge  pulses 
that  occur  in  highly  nonuniform  electric  fields  are  generally  stochastic  processes  in  the 
sense  that  when  a discharge  pulse  ceases  it  leaves  behind  an  imprint  that  can  influence 
development  of  the  subsequent  discharge  pulse.  This  imprint  may,  for  example,  be 
in  the  form  of  a negative-ion  space-charge  moving  in  the  gas,  molecules  in  excited 
metcistable  states,  or  charges  on  surfaces.  It  is  therefore  the  nature  of  the  phenomenon 
that  significant  correlations  can  exist  between  successive  discharge  pulses,  e.g.,  one 
cannot  assiune  that  pulse  amplitude  and  time  separation  from  the  previous  pulse  will 
be  independent  random  variables. 
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An  indication  of  the  stochastic  behavior  of  partial  discharge  pulses  is  not  readily  ob- 
tained from  present  commonly  used  measurement  methods.  Earher  partial-discharge 
measurement  schemes  utilizing  electrical  detection  have  been  limited  to  either  dis- 
charge pulse-height  analysis  [63-65],  pulse-frequency  or  period  measurements  [65-67], 
or  broad-band  pulse-shape  determinations  [68-70],  none  of  which  give  any  indication 
of  the  stochastic  natme  of  the  phenomenon.  A meaningful  interpretation  of  discharge 
pulse-height  distributions  often  used  for  evaluating  insulation  performance  is  impos- 
sible without  information  about  the  degree  of  correlation  among  successive  pulses. 

An  electrical  method  of  discharge  pulse  measurement  is  described  here  from  which 
direct  experimental  determinations  of  a set  of  various  conditional  pulse-height  and 
time-separation  probabihty  distributions  can  be  made.  As  will  be  shown  here,  these 
distributions  provide  a quantitative  indication  of  the  degrees  of  correlation  among 
successive  pulses  that  allow  one  to  evaluate  the  extent  to  which  the  occurrence  of 
any  discharge  pulse  affects  the  development  of  subsequent  pulses.  Examples  are 
shown  of  data  obtained  by  this  method  for  Trichel  pulses  [71,  72]  in  air,  N2/O2  and 
SF6/O2  gas  mixtures.  The  results  clearly  demonstrate  a high  degree  of  correlation 
among  successive  discharge  pulses  as  indicated  for  example  by  a strong  dependence  of 
discharge  pulse  amplitude  on  time  separation  from  the  previous  pulse.  The  observed 
stochastic  behavior  is  consistent  with  current  physical  descriptions  of  the  Trichel- 
pulse  phenomenon  [72-74],  and  illustrates  the  difficulties  that  can  be  encountered  in 
attempting  to  understand  measured  pulse-height  distributions  without  information 
about  correlation  effects  [75]. 

The  measurement  system  described  here  should  be  apphcable  to  investigations  of 
the  stochastic  behavior  in  other  types  of  pulsating  electrical  discharge  phenomena. 
Clear  evidence  of  successive  discharge  pulse  correlations  are  found,  for  example,  in  the 
operation  of  Geiger  counters  [76],  prebreakdown  current  in  dielectric  liquids  [77],  and 
for  positive  burst-type  corona  in  electronegative  gases  [64].  Partial  discharge  pulses  in 
cavities  [78]  within  sohd  insulating  materials  should  also  be  inherently  stochastic  since 
the  occurrence  of  a discharge  pulse  will  invariably  leave  behind  charged  surfaces  in  the 
cavity  which  reduce  the  local  electric  field  below  the  critical  value  required  to  sustain 
a discharge.  Some  dissipation  of  this  surface  charge  will  then  be  required  before 
another  partial  discharge  pulse  can  develop.  The  development  of  a discharge  pulse 
will  therefore  be  affected  by  residual  surface  charges  from  previous  pulses.  Similar 
argmnents  might  apply  to  pulsating  dielectric  barrier  discharges  such  as  used  for 
ozone  generation  [79]. 

A detailed  description  of  the  measurement  system  and  associated  electronic  circuitry 
is  given  elsewhere  [29]  and  wjll  not  be  covered  here.  In  this  report  we  consider  only 
the  basic  theory  of  the  measurement  concept  with  a brief  description  of  the  system. 
Examples  of  results  obtained  for  pulsating  negative  corona  discharges  (Trichel  pulses) 
in  electronegative  gas  mixtures  can  also  be  found  in  Ref’s  [30-32]. 
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Table  3.  Conditional  pulse- height  and  time-separation  distributions  up  to  order  2 


2nd 

Order 


pi{qi  j Ati-j-i)  Pi  (At,  I qi.j) 
pi{qi  I qi-j-i)  Pi  (At,  I At,_j_i) 

i>2,i>j>0 


P2{qi  I q,-j-i,qi-(-i)  P2(At,- 1 At,_j_i,  9,_fc) 

P2(At,  I At,_j_i,  At,_^_i)  p2(At,- 1 qi-j,  qi-t) 

i>2,i>j>0,i>k>0,i>t>0,t:^j 


( PoiQt)  ' 

0th  1 

. 1st  ^ 

Order  j Po{Ati) 

Order 

i>0 

>2(9. 

I 


2.3.2  Theory 

In  this  section  we  define  the  various  conditional  pulse-height  and  time-interval  distri- 
butions that  can  be  measured  by  the  method  described  here  and  indicate  how  they 
can  be  used  to  determine  the  degree  of  correlation  among  successive  discharge  pulses. 
The  basic  theory  of  conditional  probabifities  upon  which  the  following  discussion  is 
based  can  be  found  in  several  texts  (see  for  example  Ref’s  [80]  and  [81]). 

The  variables  of  particular  interest  here  are:  1)  the  discharge  pulse  magnitude,  as 
measured  by  its  amplitude  expressed  in  units  of  picocoulombs  (pC)  and  2)  the  time 
separation  between  successive  pulses,  A<„,  as  measured  in  microseconds  {ps).  It  is  as- 
sumed that  the  discharge  phenomenon  consists  of  a finite  sequence  of  m pulses  denoted 
by  the  set  of  pulse  amplitudes  and  corresponding  time  separations  {gi, Atn-i}m, 
where  n = 2,3, ...  ,m  and  m ^ 1.  The  time  separation,  Afn,  is  measured  between 
the  nth  and  (n  -H  l)st  pulses. 

The  complete  sets  of  possible  conditional  pulse-height  and  time-interval  probabihty 
distribution  functions,  up  to  second  order,  are  given  in  Table  3.  The  0th  order 
functions  are  obviously  unconditional  and  defined  such  that,  for  example,  Po{qn)dqn 
is  the  probabihty  that  any  discharge  pulse  occurs  with  amplitude  in  the  range  to 
9n  + dqn  independent  of  n.  The  distributions  pi  and  p2  are  conditional,  e.g.,  pi{qn  \ 
Atn-i)dqn  is  the  probabihty  that  the  nth  discharge  pulse  for  arbitrary  n will  have  an 
amphtude  between  and  qn  + when  its  time  separation  from  the  preceding  pulse 
is  restricted  to  the  value  A<„_i,  and  p2(9n|A<„_i,  is  the  same  with  the  added 
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restriction  that  the  amplitude  of  the  preceding  pulse  is  also  fixed  at  a value  Qn-i- 
It  should  be  noted  that  a somewhat  simplified  notation  is  used  here  for  probabihty 
distribution  functions  which  differs  from  that  often  used  in  standard  texts  (see  the 
comment  attached  to  Ref.  [80]).  It  must  also  be  emphasized  that  since  there  is  no 
way  to  select  particular  values  for  n in  the  present  measurement  system,  all  measured 
distributions  apply  only  for  arbitrary  values  of  n. 


For  the  pulsating  type  of  electrical-discharge  phenomenon  under  consideration  here, 
and  Atn  may  not  always  be  independent  random  variables.  An  immediate  indica- 
tion that  these  variables  are  not  independent  can  be  obtained  from  a discovery  that 
a measured  conditional  probabihty  distribution  is  not  equal  to  the  corresponding  un- 
conditional distribution,  e.g.,  if  it  is  found  that  Pi(gn  | Atn-i}dqn  7^  Po(gn)dgn  for  all 
allowed  Atn-i  such  that  po(At„_i)  7^  0,  then  and  Atn-i  are  dependent  random 
variables,  and  a correlation  must  exist  between  pulse  amplitude  and  time  separation 
from  the  previous  pulse.  A measure  of  the  degree  of  correlation  between  two  quan- 
tities such  as  gn  and  At„_i  can  be  found  by  computing  the  appropriate  correlation 
coefficient.  For  the  case  of  the  variables  gn  and  At„_i,  the  correlation  coefficient, 
Rq„Atn-i^  is  defined  by  [81] 


R 


1/2 


(16) 


where 

TOO 

Sqr,=  Po{gn){gu-gn)^dgn  (17) 

Jo 

= f Po{Atn-l){Atn-l  - Atn-l)^  d{Atn-l)  (18) 

Jo 

TOO  too  

‘^qn.Aln— 1 ~ / / Po(  Af„_i  )pi  (^n  j At n— 1 )(9n  1 Atn— 1 ) ^(  At^— 1 ) 

Jo  Jo 

(19) 

= / Poign)gn  dgn  (20) 

Jo 

At„_i  = / po(Atn-i)At„_i  d(At„_i).  (21) 

Jo 

Equations  (17)  and  (18)  represent  the  variances  of  gn  and  A<„_i  respectively  and 
Eq  (19)  gives  the  covariance  of  these  variables.  The  correlation  coefficient  has  the 
property  that 

1 >1  RqnAt.-.  I (22) 

where  if  Rq„,At„-i  = 0,  then  the  random  variables  gn  and  At„_i  are  uncorrelated. 

The  correlation  coefficient  can  be  computed  from  Eq’s  (16)-(21)  using  the  measured 
distributions  po(9n),  po(A<„_i),  and  pi{gn  j At„_i).  It  should  be  noted  that  Eq  (20) 
gives  the  mean  pulse  height  and  Eq  (21)  the  mean  time  interval  from  which  a repeti- 
tion rate,  l/Afn-i?  can  be  determined.  Both  of  these  quantitites  are  often  obtained 
directly  from  existing  partial-discharge  measurement  systems  [65] . 
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The  probability  that  the  nth  pulse  in  a chain  of  pulses  will  have  a particular  amplitude 
or  time  separation  from  a preceding  pulse  may  depend  in  general  on  the  amplitudes 
and  time  separations  of  all  preceding  pulses.  For  example,  the  probabihty  that  a 
sequence  of  pulses  represented  by  the  set  j = 1,2,  ...,n  — 1,  will  be 

followed  by  a pulse  with  amplitude  between  and  + dqn  is  given  by  the  following 
product  of  conditional  probabihty  functions: 

n— 1 

Yl  dqjd{Atj)  dqr,  = Po(9i)pi(Afi  1 91)^2(92  I Afi,  gi)  x 

j=i 

n — 1 

P3(Af2  I52,  Afi,5i)  • • •P2n-3(A<„_1  |9n_i,  dqjd{At j)  dq„  = (23) 

n— 1 

Po{qi)  n P2jiqj+1 1 {9.,  A<Jj)p2i-i(A<j  1 qj,  {qi,  Ati}j.i)dqjd{Atj)  dq^ 

j=i 

where,  for  example,  P2j{qj+\  \ {^i,  Af,}j)  is  the  2jth  order  conditional  distribution 
of  amplitudes  for  the  {j  + l)st  pulse  when  the  amplitudes  and  corresponding  time 
separations  for  all  previous  pulses  have  fixed  values  given  by  the  set  {5,,  At,}j. 

If  the  discharge  pulse  sequence  is  jointly  a Markov  process  [80]  for  the  variables 
and  Af„_i  in  the  sense  that  all  pulses  in  the  sequence  are  correlated  only  with  their 
immediate  predecessors,  then  the  right  hand  side  of  Eq  (24)  simplifies  to: 

n— 1 

Po(gi)/>i(Ati|gi)  n P2{^ij\qj,^tj-\)P2{qj\^tj-\,qj-\)dqjd{Atj)dqn.  (24) 

j=2 

An  indication  of  whether  or  not  an  observed  pulsating  discharge  phenomenon  behaves 
hke  a Markov  process  can  be  obtained  by  the  measurement  system  described  here 
which  allows  determination  of  the  conditional  probabihty  distributions  p\{qn  I Atn-j) 
for  n > 7 = 1, 2, 3. 

If  the  pulse  sequence  is  a Markov  process,  then  it  is  required  that 

pi(9„  I Atn-j)  = Po{qn),  (25) 

for  all  allowed  Atn-j ^ i > 2,  such  that  po{Atn-j)  ^ 0.  For  most  partial  discharge  or 
pulsating  corona- discharge  phenomena,  it  can  be  expected  that  correlations  between 
adjacent  pulses  will  overwhelmingly  predominate  so  that  to  a good  approximation 
these  phenomena  behave  like  Markov  processes.  This  is  found  to  be  the  case  for 
Trichel-pulse  discharges  in  gases  which  will  be  considered  later  as  an  example  of  an 
apphcation  of  the  method  introduced  here. 

Independent  of  whether  or  not  the  assiunption  of  a Markov  process  appHes,  it  is  nec- 
essary that  the  measured  distributions  satisfy  the  following  integral  relationships  [80] 

Po{qn)  — f Po(Af„_i)pi(5n  I Afn— 1)  ^(Afn— 1) 

Jo 


(26) 
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Po{Qn)=  / Po{qr^-l)Pl{f^tn-\\qn-.)P2{qn\Atr,-uqr^-l)dqr^-1d{Atr^-l) 
Jo  Jo 

roo 

Po{Atn)  = / Po{qn)Pi{At„  I qn)  dq^ 

Jo 

roo 

pi{qn\Atn-i)  = / po{qn-i)Pi{At„-i\qn-i)p2{qn\At„.i,q„-i)dqn^i, 

Jo 


(27) 

(28) 
(29) 


which  follow  from  the  law  of  total  probabihties.  If,  for  example,  measurements  are 
made  of  the  three  distributions  po{qn)i  Po(Af„),  and  pi{qn  \ then  it  is  desir- 

able to  perform  the  integration  indicated  in  Eq  (26)  to  show  that  it  agrees  with  the 
measured  po{qn)  and  therefore  that  the  data  on  the  measured  distributions  are  self 
consistent. 


With  the  system  described  here  it  is  also  possible  to  measure  the  conditional  proba- 
bihty  distribution  pi{qn  | A<i,2)  not  indicated  in  Table.3,  where  A<i,2  = A<„_i-f  Af„_2. 
This  distribution  satisfies  the  condition 

pi{qn  I A<i,2)  = C f Po(A<i,2  - r)pi(r  | Afi,2  - r)p2{qn  \ Afi,2  - r,  r)  dr,  (30) 
Jo 

where  pi  and  p2  in  the  integral  correspond  respectively  to  pi(A<n-i  | Atn-2)  and 
P2{qn\Atn-2iAtn-\)  obtained  from  Table  3 and  C = po(A<i,2)  is  constant  for  each 
A<i,2  which,  although  it  is  not  directly  measured  by  the  present  method,  is  of  no  rele- 
vance since  only  the  profiles  of  arbitrarily  normalized  pi{qn\  At  1,2)  are  to  be  compared 
with  integrals  on  the  right-hand  side.  If  it  can  be  assumed  that  the  phenomenon  is  a 
Markov  process,  then 

P2(9n|A<i,2  - r,r)  = pi(g„|T),  (31) 

for  all  allowed  Afi,2-  If  correlations  between  successive  time  intervals  At„_i  and 
Atn-2  are  also  found  to  be  weak  such  that 


Pi(r  |A<i,2-  r)  ~Po(t), 

then  Eq  (30)  becomes 

Pl(9„|  A<i,2)  ~ C/  Po(A<i,2  - T)po{T)pi{qn\T)dT. 

Jo 


(32) 

(33) 


A measurement  of  pi(g„  | Afi,2)  can  therefore  also  be  used  to  give  an  indication  of 
whether  or  not  the  discharge  behaves  like  a Markov  process. 

It  should  be  cautioned  that  Eqs  (26-27)  may  not  always  be  reasonable  approximations 
at  all  values  of  r.  In  the  case  of  Trichel-pulse  corona  discharges  discussed  below,  it  is 
found,  for  example,  that  there  is  always  a critical  minimum  time.  Ate,  following  any 
discharge  pulse  below  which  no  subsequent  pulse  can  occur.  Thus,  for  Afi,2  — t < Ate, 
neither  Eq  (32)  nor  Eq  (33)  can  hold.  It  is  found  [32],  however,  that  Eq  (33)  still 
gives  a good  representation  of  pi{qn  \ Ati,2)  provided  the  upper  hmit  of  integration  is 
replaced  by  At  1,2  — Ate  and  At  1,2  is  not  too  close  in  value  to  2 Ate- 
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Figure  18.  Measurement  system  giving  the  switch  configurations  needed  for  measure- 
ment of  the  various  indicated  pulse-height  and  time-interval  distributions. 

2.3.3  Measurement  System 

Figine  18  shows  a somewhat  simplified  block  diagram  of  the  measurement  system 
which  indicates  the  configurations  required  for  measurement  of  the  various  conditional 
probabihty  distributions.  The  distributions  that  can  be  measured  with  this  system 
are  listed  in  Table  4. 

The  circuits  associated  with  the  time-interval  control  logic,  pi{qn\Atn-2)  logic  and 
pulse  sorter  indicated  in  Fig.  18  are  unique  to  this  system.  Details  about  the  operation 
and  design  of  these  circuits  are  given  in  Ref.  [29] . The  time-to-amphtude  converters 
(TACi,i  = 1,2,3),  single-channel  analyzer  (SCAi,f  = 1,2),  digital-delay  generator 
(DDG),  and  multichannel  analyzer  (.MCA)  are  all  commercially  available  instruments. 
The  256-channel  MCA  is  driven  by  a microcomputer  which  enables  data  to  be  stored 
on  either  a floppy  or  hard  disk  for  subsequent  analysis. 

The  electrical- discharge  pulses  are  detected  electrically  using  a preamplifier  (Al)  con- 
nected to  an  impedance,  Z,  in  series  with  the  discharge  gap.  The  pulses  from  Al  can 
be  counted  and  observed  on  an  oscilloscope.  These  pulses  are  sent  to  an  amplifier 
A2,  the  gain  of  which  can  be  adjusted  so  that  the  range  of  observed  pulse  amplitudes 
is  within  the  range  acceptable  for  input  to  the  MCA. 

The  particular  type  of  discharge  phenomenon  used  to  test  this  measurement  concept 
is  the  well  known  [71-75]  point-to-plane  negative  (Trichel)  pulse  corona  in  an  elec- 
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Table  4.  Distributions  that  can  be  measured  with  the  present  system 


Poiqn) 

Unconditional  pulse-height  and 

Po(At„) 

4 

time-interval  distributions 

Pi  {fin  I 1 ) 

Pl{Qn  I 2) 

I A<„_3) 

Pl(9nkn-l) 

P2(9nl 


Conditional  pulse-height 
distributions. 


Pi{At„\qn) 

Pi  {Atfi  I At n— 1 ) 


Conditional  time-interval  distributions 
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tronegative  gas.  The  discharge  gap  configuration  is  indicated  in  Fig.  18  and  consists 
of  a stainless  steel  point-plane  electrode  gap.  The  point  electrode  was  maintained 
at  a constant  negative  voltage  with  respect  to  the  plane  and  the  discharge  pulse 
rate  was  controlled  by  ultra-violet  (UV)  irradiation  of  the  point.  Calibration  of  the 
pulse- amphtude  measurement  in  terms  of  charge  (pC)  was  performed  by  the  method 
previously  described  [64]. 

By  appropriate  positioning  of  the  switches  shown  in  Fig.  18,  it  is  possible  to  configure 
the  system  to  measure  the  various  indicated  distributions  directly.  Data  from  all  mea- 
surement configurations  are  acctunulated  by  the  MCA  which  contains  an  analog-to- 
digital  converter  for  pulse- height  analysis.  For  measurement  of  Poiqn),  the  discharge 
pulses  from  A2  are  fed  directly  to  the  MCA  through  gate  G3  which  is  kept  open  by 
proper  positioning  of  S4.  The  measurement  of  po{Atn)  requires  use  of  a pulse  sorter 
connected  to  two  time-to-amphtude  converters  (TACl  and  TAC2)  each  of  which  gen- 
erates an  output  pulse  of  amplitude  proportional  to  the  time  separation  between  the 
indicated  start  and  stop  input  pulses.  The  outputs  of  TACl  and  TAC2  are  fed  to  the 
MCA  through  the  OR  gate  G2  for  pulse-height  analysis.  This  arrangement  permits 
recording  of  every  sequential  time  interval,  provided  all  time  intervals  are  greater 
than  the  inherent  TAC  reset  time,  Atr,  {Atr  = 50 ps  for  the  instruments  used  here). 
The  amphtudes  of  pulses  from  the  TAC’s  are  calibrated  on  a time  scale,  and  each 
TAC  can  be  adjusted  to  cover  different  time  ranges.  For  proper  operation,  TACl  and 
TAC2  should  be  set  to  the  same  range. 

The  start  pulse  from  the  sorter  for  TACl  serves  as  the  stop  pulse  for  TAC2  and  vice 
versa.  It  is  shown  [29]  that  if  there  exist  significant  correlations  between  successive 
time  intervals  A<„  and  Afn_i,  then  it  is  possible  that  po{Atn)  will  be  distorted  if 
Atn  < Atr  for  any  n.  To  help  determine  if  this  is  a problem,  a random  sampling  of 
time  intervals  can  be  performed  using  a random  pulse  generator  to  gate  the  starting 
of  TACl  and  TAC2  [30].  If  there  are  no  errors  due  to  the  finite  TAC  reset  time,  then 
po{Atn)  should  be  the  same  for  both  normal  and  random  gated  operation. 

The  measurement  of  pi{qn  \ Afn-i)  requires  use  of  a logic-controlled  digital-delay 
generator  (DDG)  which  restricts  transfer  of  pulses  from  A2  to  the  MCA  through 
gate  G3  to  a preselected  time  window.  Details  of  the  DDG  time-interval  control  logic 
circuit  are  given  in  Ref.  [29].  Basically,  the  system  is  designed  so  that  the  DDG 
sends  a pulse  of  width  ^(Af„_i)  to  open  G3  after  a time  delay  At„_i  following  a 
discharge  pulse,  provided  no  subsequent  pulses  have  occurred  within  that  interval. 
Thus  only  pulses  that  follow  an  initial  pulse  within  a time  range  A<„_i  — 6{Atn-i)/2 
to  At„_i  + 6(At„_i)/2  are  acciunulated  by  the  MCA  for  pulse-height  analysis.  To 
minimize  errors  in  the  measurement  of  pi{qn  \ A<„_i),  it  is  desirable,  as  discussed  in 
Ref.  [29],  that  the  DDG  pulse  width  be  selected  so  that  6(A<n-i)  "C  Atn-i. 

Measurement  of  the  second  order  conditional  distribution  P2{qn\  A<„_i , 9n-i ),  requires 
use  of  a single- channel  analyzer  (SCAl)  inserted  between  A2  and  the  input  to  the 
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time-interval  control  logic.  The  SCAl  restricts  the  input  pulses  that  trigger  the 
DDG  to  have  amplitudes  that  lie  within  a selected  range  from  g„_i  — S{qn-i)  to 
Qn-i  -|-  where  again  "C  qn-i  is  a desirable  condition  to  avoid  errors. 

With  the  use  of  SCA2  and  TAC3,  the  system  can  be  configured  by  the  indicated 
positioning  of  Si,  S2,  and  S3  to  directly  measure  either  pi{Atn  | qn-i)  or  pi(At„  | 
When  measuring  the  former,  SCA2  is  used  to  select  a narrow  range  of  qn-i 
values  that  will  permit  triggering  the  start  of  TAC3,  the  output  of  which  is  fed  to  the 
MCA.  The  measurement  of  pi(A<„  | At„_i)  requires  that  SCA2  select  a narrow  range 
of  pulse  amplitudes  from  the  output  of  TACl  to  start  TAC3. 

The  determination  of  pi{qn  \ 9n-i)  can  be  performed  indirectly  by  measurement  of 
P2{qn  1 At„_i,  5„_i),  pi{Atn  l^n),  and  po{qr,-i)  and  using  the  requirement  that 

-L 

Under  some  conditions  this  distribution  can  be  measured  directly  by  configuring  the 
system  to  measure  P2{qn  \ A<„_i,9„_i)  and  then  selecting  the  DDG  delay  and  time 
window  widths  respectively  to  be  smaller  and  larger  than  the  expected  minimum 
and  maximum  pulse-time  intervals.  There  may  be  cases,  however,  where  these  ad- 
justments cannot  be  made,  such  as  when  a significant  fraction  of  the  time  intervals 
exceed  the  maximum  possible  DDG  window  setting. 

The  measurement  configuration  for  determination  of  pi{qn  | ^^1,2)  requires  a rear- 
rangement of  the  time-interval  control  logic  circuit  as  discussed  in  Ref.  [29].  The 
measurement  of  pi{qn  \ Atn-2)  involves  use  of  the  Pi{qn  \ At„_2)  logic  circuit  shown  in 
Fig.  18  together  with  configuring  the  time-interval  control  logic  so  that  it  would  other- 
wise allow  measurement  of  pi{qn  | Atn-i).  It  is  also  possible  to  measure  p\{qn  \ Atn-3) 
by  using  the  pi{qn  \ Af„_2)  logic  circuit  together  with  having  the  time-interval  control 
logic  set  for  a pi{qn  \ ^<1,2)  measurement.  In  the  case  of  a pi{qn  \ At„_2)  measurement, 
the  indicated  pi{qn  \ Atn-2)  logic  circuit,  when  inserted  by  switch  S4  in  Fig.  18,  acts  to 
inhibit  opening  of  gate  G3  until  after  there  has  been  a coincidence  between  the  digital 
delay  generator  pulse  selected  by  the  time-interval  control  logic  and  a discharge  pulse 
from  A2,  i.e.,  imtil  after  there  has  been  the  occurrence  of  a pulse  time-interval  Atn-2- 
The  next  pulse  to  occur  after  this  coincidence  is  then  allowed  to  be  recorded  and 
analyzed  by  the  MCA  through  G3  independent  of  its  time  separation,  Af„_i,  from 
the  previous  pulse.  After  a pulse  height  is  recorded,  the  MCA  sends  an  “event  pulse” 
to  reset  the  pi{qn  \ Atn-2)  logic  circuit. 

2.3.4  Results — Trichel  Pulses  in  N2/O2 

Examples  of  results  obtained  for  UV-sustained  negative  point-plane  Trichel-pulse 
corona  [71-75]  in  an  N2/20%02  gas  mixture  are  presented  here.  Figure  19  shows 


Pl{Qn\qn-l) 


P2{qn  I Af„_i , qn-\ )p\  {Atn-\  | ^n-l ) Po{qn-l  )d(  A<„_i  ).  (34) 
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data  on  the  measured  distributions  po(9n),  Pi{qn  | and  P2{qn  | for 

Trichel  pulses  (data  points)  together  with  results  of  the  integrals  given  by  Eq’s  (26) 
and  (29)  for  Po{qn)  and  pi{qn  | A<„_i)  (dashed  and  dot-dashed  lines).  These  data  were 
obtained  for  an  absolute  gas  pressure  of  100 kPa  (~  1 atm),  a point-to-plane  gap 
spacing  of  2.08  cm,  and  a point  electrode  radius- of- curvature  at  the  tip  of  0.094mm. 
Both  electrodes  were  made  of  polished  stainless  steel.  Since  the  discharge  was  UV- 
sustained,  the  pulse  rate  (related  to  mean  time  interval  Atn  from  Eq.  (21))  was 
controlled  by  the  intensity  of  the  UV  radiation  directed  at  the  point  cathode  through 
a quartz  window.  The  discharge  pulses  in  this  case  were  initiated  by  photo-electrons 
released  from  the  cathode,  and  because  the  electron  release  rate  is  roughly  propor- 
tional to  the  UV  intensity,  At„  should  decrease  with  increasing  UV  intensity.  This 
trend  is  evident  from  the  measured  data  for  po{At„)  shown  in  Fig.  20  that  indicate 
a much  broader  distribution  at  the  lower  UV  intensity  than  at  the  higher  intensity. 
The  high  UV-intensity  data  in  Fig.  20  correspond  to  the  same  conditions  under  which 
the  data  shown  in  Fig.  19  were  obtained. 

The  distributions  plotted  in  both  Fig’s  19  and  20  have  again  been  arbitrarily  normal- 
ized to  the  maxima.  Values  selected  for  At„_i  and  g„_i  used  to  obtain  the  conditional 
pulse-height  distributions  displayed  in  Fig.  19  are  indicated  in  the  legend  with  cor- 
responding plotting  symbols.  Both  po{qn)  and  pi(g„  | A<„_i)  have  been  plotted  on 
a logarithmic  scale,  and  P2{qn  \ Afn-i,9n-i)  is  shown  on  a linear  scale.  The  pulse 
heights  are  expressed  here  in  units  of  charge  (pC).  The  solid  lines  in  Fig.  19  represent 
gaussian  profiles  that  have  been  fitted  to  the  P2{qn  \ A<„_i,5„_i)  data.  The  dashed 
and  dot-dash  lines  represent  results  of  the  convolutions  from  Eq’s  (29)  and  (26)  re- 
spectively which  required  use  of  the  corresponding  data  for  po{Atn)  shown  in  Fig.  20. 
A gaussian  representation  of  pi{qn  \ Af„_i)  was  used  in  performing  the  integration 
of  Eq.  (26)).  As  evident  from  Fig.  19,  this  assumption  leads  to  errors  in  the  cal- 
culated po{qn)  at  low  Qn  values  where  the  gaussian  representation  begins  to  become 
inadequate  due  to  increasing  asymmetries  in  the  pi{qn  \ A<„_i)  profile. 

From  the  results  presented  here,  it  can  be  seen  that  there  are  significant  memory 
effects  that  influence  development  of  Trichel  pulses.  The  data  displayed  in  Fig.  19 
clearly  indicate  that  for  arbitrary  qn 

Po{qn)  ^ Pliqn  1 Atn— l)  ^ P2{qn  | Atn— li  9n— l)? 

at  allowed  values  for  Atn-i  and  ^n-i-  Thus  it  must  be  concluded  that  qn,  Atn-i,  and 
qn-i  are  not  independent  random  variables,  i.e.,  the  most  probable  amplitude  of  any 
discharge  pulse  is  influenced  by  the  time  separation  from  and  the  amplitude  of  the 
preceding  pulse.  It  is  obvious  from  the  data  on  pi  {qn  \ Atn-i)  that  the  mean  discharge- 
pulse  amplitude  increases  with  increasing  time  separation  from  the  preceding  pulse. 
This  behavior  is  expected  if  one  allows  the  possibility  that  the  negative-ion  space- 
charge  cloud  from  the  preceding  pulse  [73]  has  a shielding  effect  that  reduces  the 
electric  field  in  the  vicinity  of  the  point  electrode.  The  electric-field  strength  near 


2.3  Stochastic  Properties  of  Partial  Discharges 


45 


Q. 


c 

< 


c 


O) 

o 


Figure  19.  Normalized  unconditional  {po)  and  conditional  {pi  and  P2)  pulse-height 
distributions  for  negative,  point-plane  Trichel-pulse  discharges  in  an  N2/20%02  gas 
mixture  at  an  appHed  gap  voltage  (14)  of  9.4 kV.  The  data  points  are  from  measure- 
ments; the  dot-dashed  curve  is  from  an  integration  of  Eq  (26);  the  dashed  curves  are 
from  an  integration  of  Eq  (29);  and  the  solid  curves  are  gaussian  profiles  fit  to  the  p2 
data.  The  conditional  distributions  were  measured  at  the  indicated  “fixed”  values  for 
A<„_i  and  Qn-i-  The  maxima  in  the  distributions  are  commonly  denoted  to  occur  at 
values  9„4niix  which  are  obviously  different  for  each  distribution. 


46 


2.  GASEO US  DIELECTRICS  RESEARCH 


1.0 


0 

E 

c 

< 

Q. 

X -1.0 

c 

^~o 

O) 

o 

-•  -2.0 


-3.0 

0 100  200  300  400  500 

Atn  (^S) 


1 1 1 — 

1 1 1 1 “I 1 

A-” 

1 

« 

• 

• 

• 

• 

• 

• 

• 

• 

— r 

n □ a 

• • - 

” o 6 

“ \ 

- 0 

• 

□ 

D ° 

- 

9, 

• Low  UV  intensity,  = 9.0  kV 

■ o 8 

o High  UV  intensity,  Vg  = 9.4  kV  - 

m 

^ ^ 

1 1 1 1 1 1 

Figure  20.  Normalized  unconditional  discharge  pulse  time-interval  distributions 
measured  at  the  indicated  apphed  voltages  (14)  and  UV  intensities  for  negative, 
point-plane  Trichel  pulse  discharges  in  an  N2/20%02  gas  mixture.  The  distribution 
measured  at  9.4  kV  for  high  UV  intensity  corresponds  to  conditions  that  are  identical 
to  those  used  to  obtain  the  data  in  Fig.  19.  The  maxima  in  the  distributions  are 
commonly  denoted  to  occur  at  values  Atn,max  which  are  obviously  different  for  the 
two  distributions. 
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Figure  21.  Correlation  coefficients  for  the  variables  and  Af„_i  as  functions  of  dis- 
charge gap  voltage  (V^)  at  two  different  indicated  levels  of  UV  intensity  for  negative, 
point-plane  Trichel-pulse  discharges  in  an  N2/20%02  gas  mixture.  The  coefficients  are 
calculated  from  Eq’s  (16-21)  using  measured  distributions  such  as  shown  in  Fig’s  19 
and  20. 

the  point  will  tend  to  increase  as  the  space-charge  cloud  moves  further  away  toward 
the  anode.  Since  the  mean  discharge-pulse  size  will  increase  with  increasing  field 
strength,  relatively  larger  pulses  can  be  expected  to  develop  as  more  time  elapses 
between  pulses.  A more  extensive  analysis  of  the  data  on  pi(gn  | in  terms 

of  ion  space-charge  motion  in  the  discharge  gap  have  shown  that  the  experimental 
results  are  generally  consistent  with  this  interpretation  [32]. 

The  corresponding  correlation  coefficients,  computed  from  data  on  po(9n), 

Pi(gn  I and  po(Atn)  for  different  discharge  gap  voltages  and  UV  intensities 

are  shown  in  Fig.  21.  These  results  demonstrate  that  Rg„,At„-i  approaches  unity 
(thereby  indicating  a higher  degree  of  correlation  between  g„  and  Af„_i)  as  both  UV 
intensity  and  apphed  voltage  increase.  The  increase  in  Rg„,At„-i  with  UV  intensity  is 
reasonable  because,  as  indicated  by  the  data  in  Fig.  20,  At„  decreases  with  increasing 
UV  intensity,  thus  allowing  less  time  for  “memory  effects”  in  the  gap,  e.g.,  ion  space 
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charge,  to  dissipate.  The  increase  in  with  appHed  voltage  is  also  undoubtedly 

associated  with  an  observed  decrease  in  At„  with  voltage. 

The  data  on  P2{qn  | Atn-i,9n-i)  indicate  that  for  a given  large  pulses  are 

preferentially  followed  by  large  pulses  and  small  pulses  by  small  pulses.  The  reasons 
for  this  type  of  correlation  are  not  entirely  obvious,  but  might  be  related  to  the 
formation  of  neutral  metastable  species  during  the  discharge  pulse  that  are  more 
readily  ionized  by  electron  colhsions.  Enhancements  in  the  ionization  rates  for  the 
gas  due  to  the  presence  of  metastables  will  necessarily  lead  to  an  increase  in  discharge 
pulse  size.  Because  metastable  densities  can  be  expected  to  increase  with  increasing 
discharge  magnitude,  large  pulses  will  preferentially  give  rise  to  larger  subsequent  ^ 
pulses.  It  is  also  found  that  the  degree  of  correlation  between  and  decreases 

as  Atn-i  increases  which  is  expected  as  more  time  is  allowed  for  metastables  to  diffuse 
away  from  the  immediate  vicinity  of  the  point  electrode  where  ionization  occurs. 

Although  not  shown  here,  measurements  of  the  conditional  time-interval  distributions 
Pi(Af„  I A<„_i)  and  pi(A<„  | qn)  indicate  that  for  Trichel  pulses  there  are  measurable, 
albeit  weak,  correlations  between  A<„  and  Af„_i  and  between  Af„  and  qn-  The 
results  from  these  measured  conditional  time-interval  distributions  can  be  interpreted 
in  terms  of  the  influence  of  ion  space  charge  and  metastables  on  the  probabiHty  for 
initiating  the  next  discharge  pulse  [32]. 

It  is  also  found  from  measurements  of  pi{qn  \ Af„_2)  that  qn  and  Atn-2  are  not  signif- 
icantly correlated,  i.e.,  for  allowed  A<„_2,  Pi{qn  | Atn-2)  does  not  differ  signiflcantly 
from  poiqn)-  This  means  that  significant  memory  effects  do  not  extend  back  in  time 
beyond  the  most  recent  pulse  and  therefore  the  Trichel-pulse  phenomenon  behaves 
like  a Markov  process.  It  is,  in  fact,  expected  that  for  this  type  of  phenomenon,  the 
correlations  between  adjacent  pulses  will  dominate,  e.g.,  the  ion  space  charge  from 
the  most  recent  discharge  pulse  will  have  by  far  the  greatest  influence  on  development 
of  any  given  discharge  pulse. 


2.S.5  Results — Trichel  Pulses  in  Air 

As  expected,  the  stochastic  behavior  of  Trichel  pulses  in  air  is  found  to  be  nearly 
identical  to  that  of  Trichel  pulses  in  N2/O2  mixtures.  We  show  here  examples  of 
additional  results  that  indicate  how  Trichel  pulse  development  is  strongly  influenced 
by  memory  effects. 

Shown  in  Figinres  22  and  23  are  corresponding  sets  of  measured  data  on  po{q)  and 
Po{At)  respectively  at  the  different  indicated  discharge  point-to-plane  gap  voltages 
Va  for  Trichel  pulses  in  room  air.  The  point  electrode  was  made  of  polished  stainless 
steel  with  a tip  radius-of-curvature  of  0.01mm  and  the  point-to-plane  gap  spacing 
of  1.92  cm.  All  of  the  data  in  Figinres  22  and  23  were  obtained  with  the  same  level 
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Figure  22.  Measured  log(po(q)/po(gmax}}  versus  g at  the  different  indicated  gap  volt- 
ages Va  for  a fixed  UV  intensity.  Corresponding  data  for  po(A<)  and  pi(g„  | Atn-i) 
are  given  in  Figiures  23  and  26. 


50 


2.  GASEOUS  DIELECTRICS  RESEARCH 


At  (/us) 


Figlire  23.  Measured  po{At)/po{Atmax)  versus  At  at  the  different  indicated  gap  volt- 
ages Va  for  a fixed  UV  intensity.  Corresponding  data  for  po{q)  and  pi{qn  \ At„_i)  are 
given  in  Figures  22  and  26. 
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of  UV-irradiation.  For  ease  of  plotting,  the  distribution  functions  shown  have  again 
been  normalized  to  their  maxima,  e.g..  Figure  22  shows  Po{q) / PoiQmax)  versus  q,  where 
Qmax  is  the  value  of  q at  the  maximum  in  po{q). 

As  indicated  by  the  vertical  arrows  in  Figure  23,  there  always  exists  a minimum  time 
interval  Ate  below  which  no  subsequent  pulses  appear  after  a given  pulse.  This  min- 
imum interval  is  evidently  related  to  the  time  required  for  the  positive  and  negative 
ion  space  charge  resulting  from  a discharge  pulse  to  clear  the  region  near  the  point 
electrode  enough  so  that  the  local  electric  field  is  sufficiently  restored  to  allow  forma- 
tion of  the  next  pulse.  It  can  be  expected,  as  observed  here,  that  Ate  will  diminish 
with  apphed  gap  voltage  since  the  ions  will  be  swept  away  more  quickly  at  higher 
field  strengths. 

The  profiles  of  the  measured  po(At)  for  At  > At  max  nearly  exhibit  a poisson  character 
as  expected  if  this  part  of  the  distribution  is  determined  by  the  time  lags  associated 
with  field-enhanced  photoelectron  release  needed  to  initiate  the  next  pulse.  This 
effect  is  supported  by  the  data  shown  in  Figure  24  for  po{At)  measured  for  different 
UV  intensities  /,(/i  > I2  > I3),  where  it  is  seen  that  the  slope  of  the  log  (po(Af)) 
versus  At  curve  for  At  > At  max  decreases  with  decreasing 

The  data  on  po{At)  in  Fig’s  23  and  24  together  with  the  corresponding  po{q)  data 
in  Fig’s  22  and  25  clearly  show  a strong  correlation  between  qn  and  Atn-i-  Note  in 
Fig.  22  that  as  the  apphed  voltage  increases  the  mean  pulse  amplitude  q,  as  deter- 
mined by  Eq.  (20),  actually  decreases.  Similarly,  as  seen  in  Fig.  25,  q also  decreases 
with  increasing  UV  intensity.  This  is  obviously  the  type  of  behavior  that  would  be 
expected  if  qn  is  strongly  correlated  with  Atn-i.  As  the  mean  pulse  repetition  rate 
increases  there  should  be  a corresponding  decrease  in  q.  The  strength  of  the  corre- 
lation between  and  Atn-i  is  indicated  by  the  data  on  pi{qn  | Atn-i)  displayed  in 
Fig.  26  for  different  indicated  A<„_i  and  for  Va  = 9.2  kV. 

The  distributions  pi{qn  \ Af„_i)  tend  to  have  nearly  a gaussian  shape.  The  mean 
value  of  qn  can  decrease  significantly  with  Atn-i  as  illustrated  in  Fig.  27  which  shows 
a plot  of  the  qmax  for  pi{qn  \ At„_i)  versus  Atn-i  corresponding  to  the  results  for 
Po{q)  and  po{At)  presented  in  Figures  20  and  21.  The  amplitude  of  the  discharge 
pulse  will 'approach  a maximum  limiting  value,  qum,  as  the  time  separation  from  the 
previous  pulse  becomes  sufficiently  long  that  the  space  charge  completely  clears  the 
gap.  This  is  referred  to  as  the  “first  pulse  condition”  [68]  and  is  evident  from  the  data 
in  Fig.  27  as  a saturation  effect  and  from  the  po(9)  data  of  Fig.  25  as  a sharp  peak  or 
abrupt  drop-off  at  q = Features  associated  with  q — > qiim  are  apparent  in  the 
Po{q)  distribution  only  if  the  corresponding  po(Af)  distribution  is  sufficiently  broad 
so  that  a significant  fraction  of  the  observed  discharge  events  are  separated  from 
previous  events  by  times  greater  than  the  time  for  space-charge  transport  across  the 
gap. 


Log  (po(At)/Po(Atniax» 


52 


2.  GASEO  US  DIELECTRICS  RESEARCH 


At  ifJLS) 


Figure  24.  Measured  log(po(A<)/po(A<ma*))  versus  At  for  K = 10 kV  and  three 
different  UV  intensities,  Ii  > I2  > I3.  Corresponding  po{q)  data  appear:  iij  Fig.  25. 
The  electrode  gap  configuration  was  nearly  the  same  as  used  to  obtain  the  data  in 
Figs.  22  and  23. 
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Figure  25.  Measured  log(po(9)/pd(9max))  versus  g for  14  = 10  kV  and  three  different 
UV  intensitites,  I\  > I2  > I3.  Corresponding  po{At)  data  appear  in  Figure  24. 
The  electronde  gap  configuration  was  nearly  the  same  as  used  to  obtain  the  data  in 
Figs.  22  and  23. 
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q (pC) 


Figure  26.  Measured  pi{qn  \ ^in-i)lpi{qn,max  \ versus  q for  the  indicated 

Atn-i  values  in  microseconds  at  = 9.2  kV  together  with  the  corresponding  normal- 
ized unconditional  pulse  height  distribution  expressed  as  log(po(9)/po(9max))-  Also 
see  corresponding  data  in  Figs.  22  and  23. 


Log  (po(q)/Po(q  maxll 


Qmax  ^Or  Pi(qn|Atn-l)  (pC) 
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Atn-1  (ms) 


Figure  27.  Plot  of  qmax  versus  Af„_i  for  the  different  indicated  Va  corresponding  to 
the  data  in  Figs.  22,  23,  and  26.  Here  qmax  is  the  value  of  qn  at  which  Pi{qn  \ At„_i) 
is  at  its  maximum. 
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2.3.6  Results — Trichel  Pulses  in  SFq/02 

In  the  past  there  has  been  some  question  about  the  existence  of  Trichel  pulses  in 
SFe  [64,  82,  83].  The  present  measurements  were  performed  to  investigate  the  influ- 
ence of  SFe  on  the  stochastic  properties  of  Trichel  pulses  in  oxygen.  In  general,  it 
was  found  that  the  behavior  of  Trichel  pulses  in  SF6/O2  is  quahtatively  similar  in 
most  respects  to  that  found  for  other  electronegative  gas  mixtures  such  as  the  N2/O2 
mixture  discussed  above. 

Some  examples  of  measured  unconditional  pulse-height  and  corresponding  pulse-time- 
interval  distributions  {po{q)  and  po{At))  are  shown  in  Fig.  28.  In  this  case,  the 
discharge  was  only  partially  sustained  by  the  UV  radiation  so  that  when  the  UV  was 
turned  off,  discharge  pulses  still  occurred,  but  were,  on  average,  much  more  widely 
separated  in  time.  It  is  seen  that  in  the  case  where  the  time-interval  distribution 
becomes  quite  broad  (no  UV)  the  corresponding  po{q)  exhibits  a sharp  peak  at  as 
also  seen  from  the  data  for  air.  The  correlation  between  qn  and  At„_i  as  evident  from 
the  data  on  pi{qn  \ shown  in  Fig.  29  is  the  same  as  seen  for  air  and  N2/O2 

mixtures.  This  is  expected  since  the  physical  mechanism  of  space-charge  motion 
responsible  for  this  correlation  effect  is  essentially  the  same  for  all  electronegative 
gases. 

It  is  found  that  as  the  SFe  content  in  oxygen  increases,  the  mean  discharge  pulse 
height  decreases  and  the  mean  time  separation  between  pulses  increases  for  any  fixed 
gap  voltage,  gas  pressure,  and  electrode  configuration  [31].  Also,  as  SFe  content 
increases  the  discharge  pulses  appear  more  randomly  in  time,  i.e.,  there  is  a broad- 
ening of  the  pulse  time-interval  distribution.  It  is  evident  that  addition  of  SFe  to  O2 
has  the  effect  of  suppressing  the  growth  of  corona  discharges  due  to  the  increased 
electronegativity  of  the  gas.  In  “pure”  SFe  the  “Trichel”  pulses  that  occur  under  self- 
sustained  conditions  appear  to  be  quite  small  (below  20  pC ) and  exhibit  a very  broad 
time- interval  distribution.  Because  of  the  relatively  small  size  and  corresponding  spa- 
tial extent  of  the  Trichel  pulse  discharges  in  pure  SFe,  it  can  be  expected  that  their 
distributions  in  po{q)  and  po{At)  will  be  significantly  affected  by  details  of  electrode 
geometry  and  fluctuations  in  electron  release  rates,  gas  density,  etc.  The  present 
observations  and  interpretations  of  negative-corona  pulses  in  SFe  are  consistent  with 
previous  work  [82,  83].  Previously  observed  [64]  Trichel  pulses  at  over- voltages  for 
SFe  in  the  absence  of  UV  irradiation  seem  to  correspond  to  the  first  pulse  conditions 
[q  — ► giim)  because  of  the  relatively  low  pulse  repetition  rates,  i.e.,  long  mean  time 
intervals  between  pulses. 


Pjj(At)  / Po(Atmax)  ' Po^P 
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Figure  28.  Measured  Po{q) / po{qmax)  versus  q and  corresponding  PoiAt) / po{Atmax) 
versus  At  at  14.0  kV  in  02/10%SF6  mixture  imder  two  conditions  of  irradiation  (UV 
= on  and  UV  = off). 
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Figure  29.  Measured  pi{qn  | ^tn-i)fpi{qn,max  | versus  q for  the  indicated 

A<n-i  values  at  = 14.0  kV  together  with  the  corresponding  measured  (closed 
squares)  and  calculated  (dashed  line,  using  Eq.  (26))  unconditional  pulse-height  dis- 
tribution Po{q) / Poiqmax)  for  an  02/10%SF6  mixture. 


59 


3.  INTERFACIAL  PHENOMENA  IN 

LIQUIDS 


Task  04 

Edward  F.  Kelley 
Electricity  Division 

National  Institute  of  Standards  and  Technology 

3.1  Introduction 

The  objective  of  this  research  is  to  understand  the  mechanisms  of  electrical  failure  in 
practical  insulation  materials  used  in  power  systems  apparatus.  In  even  the  simplest 
insulation  system,  there  are  a large  number  of  factors  which  influence  the  electrical 
breakdown  voltage  of  that  system.  As  the  insulation  medium  becomes  more  com- 
plicated, such  as  in  liquid-solid  composite  systems,  the  number  of  factors  which  can 
change  the  breakdown  characteristics  increases  dramatically.  There  exists  a body  of 
tests  and  standards  for  quahfying  insulation  materials  for  suitabihty  in  power  ap- 
paratus. However,  there  often  are  foimd  to  be  wide  variations  in  the  results  from 
such  tests  because  there  are  many  factors  which  are  not  properly  controlled.  Such 
variations  are  indications  of  a lack  of  complete  understanding  of  all  the  factors  which 
contribute  to  the  failure  process.  It  is,  therefore,  the  goal  of  this  research  to  guide  the 
development  of  models  which  serve  to  quantify  the  breaJcdown  process  based  on  the 
physical  and  chemical  properties  of  the  insulation  media.  New  standards  can  then  be 
created  to  characterize  more  adequately  the  suitabihty  of  both  new  and  existing  insu- 
lation media.  The  primary  tools  for  providing  modeling  information  are  voltage  and 
current  measurements  coupled  with  optical  observations  in  the  form  of  light-emission 
measurements,  high-speed  photography,  or  electro-optical  measurements. 

To  this  end,  three  areas  of  research  are  presented  here:  1)  Partial  Discharges  as 
a Function  of  Pressure,  2)  Nanosecond  Breakdown  at  a Liquid-Solid  Interface,  and 
3)  An  Image-Preserving  Optical  Delay.  The  third  item  has  been  discussed  in  prior 
reports  without  presenting  the  details  of  the  apparatus.  The  details  are  presented  in 
this  report  essentially  as  submitted  for  the  42nd  Annual  Conference  of  the  Society 
for  Imaging  Scientists  and  Engineers,  Boston,  May  1989. 


3.2  Partial  Discharges  vs  Hydrostatic  Pressure 

The  importance  of  studying  and  understanding  partial  discharge  (PD)  phenomena 
cannot  be  overstressed,  PD  activity  exists  in  almost  all  power  systems  apparatus. 
It  is  also  anticipated  that  the  initiation  of  the  prebreakdown  process  is  the  same 
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Figure  30.  Equipment  diagram  for  observing  partial  discharges  associated  with  dc. 

for  steady-state  PD  as  for  impulse  brealcdown.  That  is,  the  prebreakdown  streamer 
initiation  for  impulse  or  for  PD’s  at  their  earhest  times  appear  to  be  identical.  This 
is  clearly  true  for  negative  polarity,  and  expected  to  be  true  for  positive  polarity. 
It  is  the  initiation  of  the  phenomena  which  is  important.  The  key  to  improving 
insulation  lies  in  the  prevention  of  the  initiation  of  the  streamer  or  PD.  Furthermore, 
the  modeling  required  to  rephcate  the  phenomena  is  likely  simpler  for  the  initiation 
than  for  the  growth  of  the  streamer. 

A collaborative  experiment  was  conducted  with  a group  from  the  University  of  Ten- 
nessee headed  by  Dr.  M.  0.  Pace.  Partial  discharges  in  hexane  were  investigated 
under  the  apphcation  of  dc  voltage  at  different  hydrostatic  pressures  from  27kPa 
to  0.3  MPa  absolute  pressure  (27kPa  is  slightly  above  the  vapor  pressure  of  hexane 
at  room  temperature).  This  represents  an  extension  of  previous  observations  which 
were  made  at  atmospheric  pressure  (0.1  MPa).  The  experiment  is  diagrammed  in 
Fig.  30.  A needle  with  a tip  radius  of  3 //m  is  placed  5 mm  from  a sphere.  A sensitive 
amphfier  monitors  the  current  supplied  to  the  PD,  and  an  image- converter  camera 
photographs  the  growth  of  the  PD.  In  order  to  capture  the  initiation  of  the  PD,  a 
random  event,  it  is  necessary  to  employ  an  image-preserving  optical  delay  to  store 
the  image  information  in  air  imtil  the  camera  can  be  triggered  by  the  random  event 


3.2  Partial  Discharges  vs  Hydrostatic  Pressure 


61 


Figure  31.  Partial  discharge  photographs,  from  left  to  right:  27kPa,  2fs  between 
frames;  0.1  MPa,  200ns  between  frames;  0.2  MPa,  200ns  between  frames. 


and  begin  taking  pictures. 

The  reason  for  using  hexane  is  that  a large  body  of  experimental  results  already  exists 
in  the  literature  and  the  liquid  is  well-characterized  both  physically  and  chemically. 
Liquids  like  transformer  oil  are  of  more  practical  interest  but  suffer  from  not  being 
well-charaeterized  chemically  or  physically  and  large  variations  can  be  seen  from 
barrel  to  barrel.  It  is  anticipated  that  the  results  established  for  hexane  or  any  other 
simple  liquid  will  be  readily  apphcable  to  all  liquids.  However,  at  this  early  stage  we 
want  to  avoid  as  many  complications  as  possible  which  arise  from  the  selected  liquid. 

The  results  of  this  experimental  study  confirm  that  a PD  is  a gaseous  chamber  - a 
bubble.  In  general,  the  size  of  the  bubble  decreased  with  an  increase  in  pressure, 
and  the  amoimt  of  charge  injected  during  the  expansion  of  the  bubble  decreased  as 
the  pressure  increased.  Figure  31  shows  PD  growth  at  several  pressures.  The  rings 
observed  in  each  frame  have  nothing  to  do  with  the  experiment;  they  arise  from  the 
interference  of  the  laser  beam  at  the  surface  of  the  image-converter  camera  tube.  The 
region  viewed  by  the  camera  in  the  direction  the  needle  is  pointing  is  0.29  mm  from 
one  edge  of  the  frame  to  the  other.  The  smallest  particle  which  could  be  discerned 
in  the  photograph  is  estimated  to  be  5 fim  in  diameter.  The  growth  of  the  PD  near 
the  vapor  pressure  of  hexane  is  continuous,  that  is,  the  PD  bubble  keeps  growing, 
eventually  breaking  up  into  several  smaller  bubbles  which  leave  the  area  of  the  gap 
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and  are  readily  visible  to  the  eye.  For  higher  pressures,  the  PD  is  seen  to  collapse 
back  on  itself  with  at  most  a few  very  small  bubbles  leaving  the  region  of  the  gap. 
These  small  bubbles  leave  the  region  of  the  gap  and  quickly  dissolve  into  the  liquid. 

The  initiation  of  PD’s  under  the  apphcation  of  impulses  has  been  studied  before  in 
nitrobenzene  [84].  A hypothetical  model  for  streamer  initiation  was  proposed:  There 
is  a favorable  site  on  the  tip  of  a needle  at  which  charge  injection  is  large.  As  the 
field  increases,  the  “jet”  of  charge  injection  increases  until  the  joule  heating  of  the 
liquid  causes  a bubble  to  initiate.  The  bubble  grows  and  branches  becoming  the 
streamer.  The  rapid  changing  of  a small  portion  of  a liquid  into  a gas  represents 
a small  explosion  in  the  liquid.  The  initial  bubble  would  be  at  high  pressure.  This 
hypothesis  explains  the  observation  of  an  initial,  spherical  shock  wave  emanating  from 
the  tip  of  the  needle  at  the  moment  of  initiation  as  observed  in  many  photographs  of 
streamer  growth  under  the  apphcation  of  impulses.  To  what  extent  this  picture  can 
be  apphed  to  PD’s  under  the  apphcation  of  steady-state  voltages  remains  to  be  seen. 

Attempts  have  been  made  to  apply  various  models  to  the  growth  of  a PD  [85].  Wat- 
son [86],  Alexeff  [87],  and  Fenimore  [88]  have  developed  models  which  are  based  on 
Coulombic  repulsion  of  the  charges  on  the  surface  of  the  bubble  or  the  heating  of  the 
gas  within  the  bubble  due  to  the  discharges  as  the  bubble  expands.  An  additional 
complication  would  be  the  initial  explosion  of  the  liquid  into  a gas  at  the  initiation 
site.  It  is  difficult  to  determine  exactly  what  process  is  causing  the  expansion  of 
the  PD  bubble.  It  may  be  a combination  of  the  initial  explosion,  Coulombic  repul- 
sion, heating  of  the  gas  within  the  bubble,  and  further  vaporization  of  liquid  at  the 
boundary  of  the  bubble.  Each  of  these  processes  may  play  a role,  and  the  respective 
importance  of  each  process  may  change  as  the  PD  develops  into  a streamer  and  the 
streamer  grows  across  the  gap.  Only  further  quantitative  results  will  provide  insight 
to  distinguish  or  combine  these  pictures. 

Regarding  future  developments:  At  the  beginning  of  the  experiment,  attention  was 
paid  to  particulate  content  of  the  hexane.  Ultrapure  hexane  was  used  at  first  in 
combination  with  circulating  the  hexane  through  a 1 — fim  polytetrafluoroethylene 
particle  filter.  With  such  purity,  it  was  difficult  to  observe  any  partial  discharges. 
When  the  cell  system  was  refilled  with  spectroscopic  grade  hexane  and  no  particle 
filtration  was  attempted,  partial  discharges  were  readily  observed.  This  suggests 
that  small  particles  play  an  important  role  in  the  generation  of  partial  discharges. 
To  explore  this  area  further,  a redesigned  high-magnification  optical  system  will  be 
needed  in  order  to  provide  micron  or  sub-micron  optical  resolution.  Such  a system 
will  be  implemented  next  year.  In  addition,  a calibration  of  the  amplifier  system  will 
be  made  to  determine  the  charge  deposited  in  the  PD  during  its  growth.  It  is  hoped 
that  these  data  will  provide  some  quantitative  measure  to  guide  the  development  of 
a model  for  PD  generation  and  growth. 
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BREAKDOWN  VOLTaGE  VS  TIME  TO  CROSS  GAP 


Figure  32.  Liquid  breakdown  data  showing  the  decrease  in  breakdown  voltage  as  the 
time-to-breakdown  is  reduced. 

3.3  Nanosecond  Breakdown  At  A Liquid- Solid  Interface 

In  previous  reports  it  has  been  shown  that  breakdown  characteristics  in  the  sub- 
microsecond  regime  differ  from  those  in  the  microsecond  regime.  Specifically,  the 
breakdown  voltage  does  not  increase  steadily  as  the  time-to-breakdown  decreases.  In 
the  case  of  a needle-sphere  electrode  system  with  a 5 mm  gap  in  transformer  oil,  as 
the  breakdown  time  decreases  to  the  range  of  a few  microseconds,  the  breakdown 
voltage  begins  to  decrease.  When  the  breakdown  times  become  significantly  shorter 
than  one  microsecond,  the  breakdown  voltage  again  increases  as  the  breakdown  time 
decreases,  as  seen  in  Fig.  32.  An  explanation  was  proposed:  As  the  voltage  increases, 
the  streamer  spends  less  time  in  the  slower  propagation  modes.  If  the  voltage  increases 
fast  enough,  the  slower  modes  may  be  circumvented  almost  entirely  and  the  streamer 
is  exclusively  comprised  of  fast  phenomena.  It  is  because  of  this  faster  mode  that  the 
breakdown  voltage  decreases  slightly  as  the  time-to-breakdown  is  made  to  decrease. 
However,  the  decrease  is  short-lived,  and  the  breakdown  voltage  again  increases  dra- 
matically when  the  streamer  is  entirely  in  the  fast  mode  and  the  time-to-breakdown 
is  forced  to  decrease  further. 

It  is  believed  that  the  faster  propagation  modes  are  more  energetic  in  some  sense 
than  the  slower  modes.  A piece  of  paper  insulation  may  stop  the  slower  mode,  but 
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Figiire  33.  Experimental  configuration  of  paper  interface  placed  in  the  gap  between 
the  needle  and  the  sphere  electrodes. 


it  may  not  stop  the  faster  mode.  If  the  fast  mode  can  be  initiated  early  in  the  life  of 
the  streamer  due  to  the  apphcation  of  a fast,  nanosecond  rise-time  pulse,  the  paper 
insulation  may  fail  at  a lower  voltage  than  for  a slower,  microsecond  pulse  where  only 
a slow  process  would  be  generated.  It  is  the  purpose  of  this  experiment  to  shed  some 
light  on  the  difference  of  fast  and  slow  streamer  processes  in  their  interaction  with 
interfaces. 

The  present  experiment  adds  a paper  interface  between  the  needle  electrode  and  the 
spherical  electrode  (a  strip  13  mm  wide  having  a thickness  of  64/[fm).  The  position 
of  the  paper  in  the  gap  can  be  varied  relative  to  the  needle  and  the  sphere,  (see 
Fig.  33).  An  image- converter  camera  is  used  to  photograph  the  time-evolution  of 
the  streamers  in  the  gap,  and  an  image-preserving  optical  delay  can  be  employed  to 
temporally  isolate  the  faster  phenomena.  We  were  unable  to  determine  whether  or 
not  a fast  pulse  can  break  down  an  interface  at  a lower  voltage  than  a slow  pulse 
can.  There  were  two  diflBculties:  1)  The  slower  streamers  existed  long  enough  that 
they  could  spread  to  the  outer  edge  of  the  paper  and  then  jump  across  the  gap,  and 
2)  the  breakdown  voltages  needed  to  completely  test  the  hypothesis  exceeded  the 
300 kV  peak  capabilities  of  the  present  system,  (see  Fig.  34).  Thus,  more  work  will 
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WITHSTAND  VOLTAGE  WITH  AND  WITHOUT  INTERFACE 


Figure  34.  Increase  in  brealcdown  voltage  with  interface  present  using  7 fis  wide  pulse. 


be  necessary  to  determine  the  answer.  However,  there  is  some  evidence  that  fast 
pulses  can  produce  anomalously  low  failures  of  an  interface. 

Figure  35  shows  three  photographic  sequences.  The  left  two  are  taken  under  identical 
voltage  conditions  of  a 7 fzs  width  pulse  having  a peak  of  103  kV.  The  only  difference 
is  the  position  of  the  interface.  The  left  photograph  shows  a fast  streamer  quickly 
penetrating  the  interface  which  is  far  away  from  the  tip.  In  the  middle  photograph, 
the  slower  streamer  process  is  stopped  by  the  paper  interface  which  is  near  the  tip  of 
the  needle.  The  right  photograph  shows  the  penetration  of  the  paper  by  a fast  event 
where  the  camera  operated  at  ten  times  the  hraming  rate  of  the  other  two  photographs 
and  the  pulse  had  a peak  voltage  of  178kV.  In  the  right  photograph  the  liquid  was 
purposely  made  dirty  by  generating  many  breakdowns  in  the  oil  thus  producing  a 
large  number  of  fine  carbon  particles  suspended  in  the  oil.  The  difference  between 
the  fast  event  in  the  left  photograph  and  the  fast  even  in  the  right  should  be  noted: 
the  particles  tend  to  reduce  the  number  of  branches  of  the  fast  event  and  may  speed 
it  up.  The  extent  of  the  slow  event  varies  from  shot  to  shot  before  the  onset  of  the 
fast  event.  In  general,  the  time  spent  in  the  slower  mode  is  reduced  with  the  addition 
of  particles  or  increasing  the  rise-time  of  the  voltage  pulse. 

What  is  clear  from  these  results  is  this:  1)  If  a fast  streamer  is  generated,  it  readily 
punctures  the  paper  insulation.  2)  The  presence  of  dirt  enhances  the  probabihty  of 
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Figure  35.  Photographic  sequences  showing  interaction  of  streamers  with  the  paper 
interface.  Two  left  photographs  have  a frame  spacing  of  500  ns,  the  right  photograph 
has  a frame  spacing  of  50ns. 
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initiating  a fast  streamer  and  the  speed  of  the  fast  streamer  is  increased.  3)  Fast 
streamers  are  very  difficult  to  initiate  in  clean  liquids  without  first  producing  the 
slow  streamer  process.  How  any  of  these  factors  combine  in  a practical  insulation 
system  remains  to  be  seen.  If  partial  discharge  activity  is  present  in  power  apparatus, 
then  there  are  particles  being  generated  and  the  above  results  would  indicate  an 
increased  susceptibility  to  anomalous  failure  due  to  fast  pulses  which  generate  the 
fast  streamer  processes.  However,  more  work  is  needed  to  isolate  the  parameters 
which  may  combine  to  permit  an  anomalous  failure. 


3.4  An  Image- Preserving  Optical  Delay 

An  image-preserving  optical  delay  (IPOD)  is  described  which  delays  an  image  from 
reaching  a high-speed  camera  for  several  hundred  nanoseconds.  This  device  permits 
the  photography  of  phenomena  preceding  fast,  random  events  by  using  the  event  as 
a trigger  for  a high-speed  camera. 

The  need  for  delayed  images  arises  when  information  is  desired  concerning  the  de- 
velopment of  very  fast,  random  phenomena  having  durations  from  nanoseconds  to 
microseconds.  Consider  a random  event  which  occurs  during  an  interval  of  time.  The 
longer  the  time  interval,  the  more  difficult  it  is  to  synchronize  a high-speed  camera 
with  the  event  under  study.  Delaying  the  image  from  reaching  the  camera  permits 
the  camera  to  photograph  events  back  in  time.  If  the  image  delay  is  long  enough, 
the  camera  photographs  phenomena  which  preceded  the  event  triggering  the  camera. 
Several  ways  have  been  successfully  employed  to  delay  an  image  [89,  90].  This  paper 
presents  another  alternative  which  has  the  potential  of  producing  delays  up  to  a few 
microseconds  in  a relatively  small  space. 

Details  of  the  Apparatus: 

Figure  36  shows  the  use  of  the  IPOD  in  an  experiment  to  observe  photographically 
electrical  brealcdown  in  liquids.  Such  breakdown  processes  are  examples  of  fast,  ran- 
dom phenomena  like  lightning.  A high-voltage  (HV)  pulse  having  a duration  of  7 fs 
is  apphed  to  a needle-sphere  electrode  system  with  a gap  separation  of  7 mm.  The 
electrodes  are  placed  within  a cell  (not  shown)  filled  with  liquid  hexane.  A xenon 
flashtube  provides  illumination  to  permit  shadowgraph  photography  using  an  image- 
converter  camera  which  records  the  time  evolution  of  the  brealcdown  process.  At  some 
random  time  during  the  application  of  the  electrical  pulse,  a streamer  will  begin  at  the 
tip  of  the  needle  and  grow  across  the  electrode  gap.  When  the  streamer  reaches  the 
grounded  spherical  electrode  breakdown  occurs.  The  camera  is  triggered  by  means  of 
a pulse  transformer  which  senses  the  high  current  pulse  of  several  hundred  amperes  in 
the  ground  line  generated  by  the  breakdown.  While  the  camera  is  preparing  to  take 
the  first  frame,  the  image  of  the  process  occurring  in  the  electrode  gap  is  delayed  by 
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Figure  36.  Image- preserving  optical  delay  orientation  in  a liquid  breakdown  experi- 
ment. 


the  IPOD  so  that  when  the  camera  begins  taking  pictures,  it  records  prebreakdown 
images. 

The  IPOD  is  the  arrangement  of  mirrors  between  the  object  (the  region  between  the 
electrodes)  and  the  camera.  Its  main  elements  are  an  array  of  flat  mirrors  and  a large- 
diameter  concave  spherical  mirror.  The  image  of  the  object  is  repeatedly  refocused 
and  colhmated  by  the  interplay  of  the  flat  mirror  array  and  the  spherical  mirror.  The 
operation  of  the  IPOD  is  best  understood  by  considering  an  equivalent  lens  system 
each  of  which  has  a focal  length  f.  When  an  object  is  placed  a distance  2f  from  a lens, 
an  unmagnified  image  will  be  produced  at  a distance  2f  on  the  other  side  of  the  lens. 
If  another  lens  is  placed  at  the  location  of  the  image,  the  second  lens  only  serves  to 
colhmate  the  light  from  the  first  lens.  Repeating  this  two-lens  arrangement  with  a 
separation  of  2f  provides  a means  of  relaying  an  image  over  a long  distance  with  imit 
magnification.  In  the  IPOD,  the  lenses  are  replaced  with  the  concave  mirror,  and  the 
flat  mirror  array  serves  to  direct  the  repeated  reflections. 

The  alignment  of  the  IPOD  is  accomplished  by  placing  the  center  of  the  flat  mirror 
array  at  the  on-axis  focal  point  of  the  concave  mirror.  Alignment  is  accomplished 
using  a He-Ne  laser  by  passing  it  through  the  electrode  gap  to  the  two  flat  entrance 
mirrors  which  steer  the  beam  into  the  array.  The  two  flat  entrance  mirrors  are 
adjusted  so  that  the  distance  from  the  object  to  the  first  mirror  of  the  array  (see 
Fig.  37)  is  equal  to  one  focal  length  of  the  concave  mirror.  Since  off-axis  imaging 
from  a mirror  introduces  astigmatism,  it  is  necessary  to  exploit  the  entire  surface  of 
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Figure  37.  Flat  mirror  array  as  viewed  from  the  concave  mirror.  The  solid  lines 
represent  the  long  path,  and  the  dashed  lines,  the  short  path.  The  central  mirror  of 
the  array  is  depicted  by  a dashed  circle. 

the  concave  mirror  in  a synunetrical  manner  so  that  the  net  astigmatism  associated 
with  the  final  image  after  the  array  is  almost  zero.  This  is  accomplished  in  the 
following  manner:  Mirror  #1  is  adjusted  so  that  the  reflected  beam  hits  the  concave 
mirror  along  a line  which  parallels  the  axis  of  the  concave  mirror.  The  return  beam 
will  then  hit  the  central  mirror  of  the  array.  The  beam  returns  to  the  concave  mirror 
and  is  reflected  back  to  mirror  ^2.  This  is  referred  to  as  a long  path.  Mirror  ^2  is 
adjusted  to  reflect  the  beam  back  to  the  concave  mirror  so  that  the  return  beam  hits 
mirror  ^3,  a short  path.  The  long  path  hits  the  concave  mirror  twice,  the  short  path 
hits  it  once.  This  sequence  is  repeated  until  the  output  mirror  #11  is  reached,  then 
the  beam  is  directed  into  the  image-converter  camera  by  the  flat  exit  mirrors.  The 
distance  between  mirror  #11  and  the  photocathode  of  the  camera  is  one  focal  length 
of  the  concave  mirror.  It  must  be  noted  that  the  short  path  to  mirror  #3  skips  the 
mirror  adjacent  to  mirror  #1.  Only  by  skipping  the  adjacent  mirror  in  this  manner 
for  each  short  path  can  it  be  assinred  that  the  beam  symmetrically  passes  across  the 
concave  mirror  equally  in  all  directions  whereby  astigmatism  is  virtually  eliminated. 

In  like  manner,  an  IPOD  can  be  made  with  6, 10, 14, 18, 22,  ...4n-f2  = N(n  = 1, 2, 3, ...) 
mirrors  in  the  circle  aroimd  the  central  mirror  and  with  one  extra  mirror  (N  -|-  1) 
outside  the  circular  array.  This  produces  3n  -|- 1 images  on  the  concave  mirror  and 
one  image,  the  final  image,  outside  the  array.  A total  of  3n  -|-  2 images  are  produced 
over  a path  length  of  (12n-f8)f.  For  such  an  array,  the  long-path-short-path  sequence 
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Figure  38.  Results  of  the  image-preserving  optical  delay,  50ns  between  frames. 

must  always  be  maintained  where  the  adjacent  mirror  is  skipped  as  described  above. 
In  the  case  of  the  IPOD  presented  here,  the  focal  length  of  the  30-cm  diameter 
concave  mirror  is  f = 372  cm;  the  flat  mirror  array  contains  10  (with  n = 2),  2.5- 
cm  diameter,  small  mirrors  around  a 5-cm  diameter  center  mirror.  The  total  path 
length  from  object  to  final  image  is  119m  which  provides  for  an  image  delay  of  397 ns. 
The  central  array  mirror  is  larger  than  the  circumferential  mirrors  to  expedite  the 
alignment  process. 

The  IPOD  functions  as  a long  focal-length  mirror  with  a small  aperture.  The  resolu- 
tion of  the  device  is  limited  not  by  the  large-diameter  mirror,  but  by  the  circumfer- 
ential array  mirrors  acting  as  apertures.  With  the  2. 5-cm  mirrors,  a resolving  power 
of  the  IPOD  can  be  calculated  and  measured  to  be  approximately  5 line-pairs  per 
millimeter  (Ip/mm).  To  improve  the  net  resolution  of  the  IPOD  without  replacing 
the  small  array  mirrors  with  larger  ones,  it  is  a simple  matter  to  magnify  the  image 
before  the  IPOD,  then  demagnify  it  by  a factor  of  two  or  more  after  the  IPOD.  The 
net  resolution  will  be  lOlp/mm  or  better  at  the  camera. 

Figure  38  shows  the  results  of  using  a premagnification  of  5.4X  before  the  IPOD  and 
a demagnification  of  2X  after  the  IPOD.  Most  of  the  frames  shown  in  the  photograph 
were  stored  in  the  air  by  the  IPOD:  The  camera  was  triggered  a few  nanoseconds 
after  the  breakdown  occurred,  at  this  point  the  phenomena  in  the  gap  appeared  as 
shown  in  frame  five.  The  camera  began  taking  pictures  approximately  150  ns  later, 
when  the  process  had  advanced  to  the  stage  shown  in  frame  eight  or  nine.  Because  of 
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the  IPOD  delaying  the  arrival  of  the  image  to  the  camera,  it  is  able  to  photograph  the 
phenomena  before  the  camera  was  triggered,  that  is,  the  IPOD  permits  the  camera 
“see”  back  in  time  several  hundred  nanoseconds. 

The  success  of  the  IPOD  depends  upon  the  quahty  of  the  mirrors  used.  Because  of  the 
many  reflections  involved,  it  is  necessary  to  use  high-reflectivity,  dielectric  mirrors. 
There  is  a reflection  at  each  distance  f for  a total  of  12n  -f  7 reflections  (n  = 2 for 
10  circumferential  array  mirrors,  see  above).  Adding  an  additional  three  entrance 
mirrors  and  three  exit  mirrors  results  in  a total  of  12n-|- 13  reflections.  If  each  mirror 
had  a reflectivity  of  97%  (accounting  for  dust,  scratches,  etc.),  the  transmittance  of 
the  IPOD  would  be  32%. 

Adding  more  mirrors  to  the  circumferential  array  will  increase  the  delay  but  reduce 
the  transmittance.  If  the  10  small  mirrors  were  replaced  with  14  (n  = 3)  mirrors  of 
5-cm  diameter  using  the  same  372-cm  focal  length  concave  mirror,  a delay  of  546  ns 
would  be  provided  with  a resolution  of  10  Ip/mm,  a resolution  commensurate  with 
or  better  than  many  high-speed  recording  devices.  However,  there  would  be  up  to 
49  reflections  (using  six  insertion  mirrors)  which  would  provide  for  a transmittance  of 
22%  through  the  IPOD  assuming  97%  reflectance  of  the  mirrors.  It  may  be  possible 
to  increase  the  delay  further  by  making  two  concentric  circumferential  arrays.  For 
example,  a 10-mirror  inner  array  with  a 16  mirror  outer  array  would  provide  for  a 
993  ns  delay.  However,  a small  amoimt  of  astigmatism  will  be  introduced  since  a way 
has  not  yet  been  discovered  to  exploit  such  a compound  array  in  a fully  symmetrical 
manner.  Such  an  array  (n  = 6,  26  circumferential  mirrors)  would  have  85  reflections 
giving  a transmittance  of  7.5%  with  a mirror  reflectance  of  97%.  The  delay  can  also 
be  lengthened  by  increasing  the  focal  length  of  the  concave  mirror.  However,  to  avoid 
a loss  of  resolution,  the  diameter  of  the  mirrors  must  be  increased  proportionately. 
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4.1  Introduction 


The  design  of  free-standing  voltage  dividers  that  are  used  to  measure  lightning  and 
switching  surges  in  power  equipment  testing  generally  result  in  dividers  of  large  phys- 
ical dimensions  in  order  to  withstand  voltages  from  several  hundred  kilovolts  to  mega- 
volts. These  devices  typically  exhibit  resonances  due  to  stray  capacitances  and  in- 
ductances associated  with  the  leads  connecting  them  to  the  voltage  source.  Problems 
arise  when  there  is  significant  energy  in  the  pulses  to  excite  these  inherent  resonances, 
distorting  the  measurement.  The  lEC  [91]  and  IEEE  [92]  standards  concerning  the 
use  of  these  dividers  for  high  voltage  impulse  measurements  recommend  evaluation 
of  their  measurement  fidelity  through  step-response  tests  and  alternatively,  through 
simultaneous  comparison  measurements  made  with  a well- characterized  reference  de- 
vice. There  is  no  guidance  provided  by  the  standards  with  regard  to  interactions 
between  dividers,  i.e.,  how  resonances  due  to  the  presence  of  one  of  the  dividers  in 
the  circuit  affects  the  measurements  made  with  the  other.  These  interactions  affect 
the  measurements  of  peak  voltage  and  consequently  the  evaluation  of  the  scale  factor 
for  the  divider  under  comparison. 

The  work  reported  here  addresses  these  problems  in  the  measurement  of  fast  voltage 
transients  through  an  assessment  of  the  magnitude  of  divider  interactions  in  simulta- 
neous measurement  configurations  and  makes  recommendations  for  their  minimiza- 
tion. This  study  is  reported  in  greater  detail  in  reference  [93].  The  results  show 
that  interaction  effects  are  most  significant  in  step-response  measurements  and  less 
important  in  the  measurement  of  standard-lightning  impulses.  The  following  sec- 
tions describe  the  purpose,  experimental  design,  and  results  of  the  study  on  divider 
interactions  in  comparison  measurements. 
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4.2  Interactions  Between  Two  Dividers  Used  in  Compari- 
son Measurements 

4.2.1  Introduction 

The  accuracy  of  an  impulse  measurement  is  determined  by  the  dynamic  behavior 
of  the  complete  electrical  system  including  the  divider,  the  test  object,  the  impulse 
generator,  interconnecting  leads,  recording  instrument,  and  chopping  gap,  if  used. 
lEC  [91]  and  IEEE  [92]  standards  concerning  impulse  voltage  measurement  tech- 
niques recommend  that  step-response  parameters  be  evaluated  as  one  approach  to 
characterize  the  measuring  system’s  dynamic  behavior.  This  approach  has  difficul- 
ties since  systems  having  step-response  parameters  within  the  limits  prescribed  by  the 
standard  can  still  have  unacceptable  errors  in  the  actual  impulse  measurements  [94] . 
Alternatively,  it  is  recommended  that  evaluation  be  made  through  a simultaneous 
comparison  measurement  between  a standard  reference  divider  and  the  one  to  be 
evaluated.  lEC  publication  60-4  [95]  calls  for  the  independent  reference  system  to 
have  overall  measurement  errors  of  less  than  1%  over  its  range  of  use.  This  study 
was  undertaken  because,  although  the  standards  specify  less  than  1%  measurement 
uncertainty,  no  indication  is  provided  concerning  the  types  of  interactions  between 
dividers  and  no  guidance  given  on  what  precautions  are  to  be  taken  to  insure  that 
they  do  not  affect  the  divider  comparison. 

Previous  comparisons  between  dividers  where  one  divider  was  used  as  a standard  to 
determine  the  accuracy  of  the  other  have  been  published  [96-99]  including  some  which 
discussed  divider  interactions.  In  this  study,  interactions  arising  from  the  simultane- 
ous use  of  two  dividers  and  resultant  distortion  of  the  signals  have  been  measured. 
Two  resistive  dividers  considered  to  be  of  reference  quahty,  a pure  capacitive  divider, 
and  a damped  capacitive  divider,  were  used.  The  comparisons  were  made  using 
both  independent  and  simultaneous  measurements  of  the  step  responses,  standard- 
lightning impulses,  and  chopped-lightning  impulses.  The  experimental  approach  is 
described  in  the  next  section. 


4.2.2  Experimental  Approach 

The  two  resistive  divider  measuring  systems  used  in  this  study  are  designated  NBSl 
and  NBS3,  rated  at  1 MV  and  0.3  MV,  respectively,  for  chopped-lightning  impulses. 
Both  dividers  have  high-voltage  arms  consisting  of  10  kf)  and  low-inductance  resistors. 
Each  has  a top  ring  to  grade  the  distributed  capacitances.  The  capacitances  for  the 
top  rings  measured  to  ground  are  69  pF  for  NBSl  and  22  pF  for  NBS3. 

Two  laboratory  model  capacitive  dividers  were  constructed  to  obtain  responses  that 
were  very  different  from  those  of  the  resistive  dividers,  NBSl  and  NBS3.  One  is  a pure 
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capacitive  divider  having  a slow  rise  time  and  46%  overshoot.  The  second  is  a damped 
capacitive  divider  with  a high-impedance  arm  consisting  of  a total  capaeitance  of 
430pF  and  total  resistance  of  5100  in  series.  Its  step-response  rises  very  slowly  and 
has  no  overshoot. 

The  experiments  utilized  two  different  pulse  generators:  a triggered-gap  generator 
that  produced  pulses  of  several  kilovolts  in  amplitude  and  one  using  silicon-controlled 
rectifiers  (SCR’s)  for  pulses  up  to  400  volts.  The  generators  produced  both  full 
standard  lightning  and  chopped-lightning  impulses.  Two  identical  scan-converter 
digitizers  with  a 400  MHz  bandwidth  and  stated  digital  resolution  of  9 bits  were  used 
for  the  step-response  measurements.  For  the  impulse  voltage  measurements,  a third 
digitizer  having  a 125  MHz  bandwidth  was  added.  Computer  programs  were  written 
to  calculate  the  time  parameters  from  the  step-response  and  lightning-impulse  digital 
data  records  according  to  the  standards  [91,  92].  The  data  were  smoothed  using  cubic 
splines  under  least-squares  constraints. 

A mercury- wetted  relay  having  a nanosecond  closure  time  was  used  with  a dc- voltage 
supply  to  generate  the  step-response  signals.  The  measuring  system  was  charged  to 
a preset  dc  voltage  and  then  discharged  to  ground  through  the  relay. 


4-2.3  Step  Response  Measurements 

Step  responses  of  the  measurement  systems  were  determined  using  the  horizontal- 
loop  and  square-loop  configurations.  For  the  horizontal-loop  configuration,  shown 
in  Figure  39,  the  signal  source  was  mounted  on  a platform  at  a height  of  2 m so 
that  the  connecting  lead  between  the  source  and  the  divider  ran  perpendicularly  to 
the  divider’s  vertical  axis.  The  source  was  connected  to  the  groimd  plane  by  a low- 
inductance  copper  strip,  labeled  G in  the  figure.  The  angular  separation  between 
the  two  horizontal  divider  leads  (^)  was  about  60  degrees  and  the  distances,  between 
dividers  {£c  and  £»)  were  3.2  m.  Preliminary  measurements  using  angular  separa- 
tions of  60,  90,  and  180  degrees  showed  that  changing  angular  separation  made  little 
difference  in  the  results. 

For  the  square-loop  configuration,  the  signal  source  was  placed  on  the  ground  plane 
and  the  connecting  lead  was  run  vertically  up  from  the  generator  to  a height  of  2 m 
and  then  across  horizontally  to  the  divider.  The  same  circuit  geometry  as  used  for  the 
horizontal  loop  configuration  was  preserved,  but  with  the  signal  source  at  a different 
location  in  the  loop.  The  horizontal  connecting  leads  were  4.6  m long  for  the  NBSl 
divider  and  3.8  m long  for  the  others.  A 50fl  attenuator  was  connected  to  measure 
the  output  of  the  signal  source  directly.  The  same  circuit  geometry  as  shown  in 
Figure  39  was  used  for  the  full-  and  chopped-lightning  impulse  measurements  but 
with  a different  signal  source. 
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Figure  39.  Experimental  horizontal  loop  configuration  used  to  investigate  simultane- 
ous measurements  using  different  dividers.  Shown  are  the  two  divider  systems  and  the 
damping  resistors  (Ra  and  Rc)  used  with  each  individual  divider.  For  the  square-loop 
configuration,  the  signal  source  is  placed  on  the  groxmd  plane. 
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Figure  40.  Step-response  of  the  independent  NBS3  divider  for  different  damping 
resistors,  (1)  Rds  = 00;  (2)  Rr^3  = 1500;  (3)  = 3500;  (4)  Ras  = 5500. 


Thirteen  different  test  configurations  were  investigated.  During  preliminary  tests, 
measurements  were  made  with  the  dividers  separated  by  1.65,  3.8,  and  8.6  m.  No 
significant  changes  in  interaction  effects  were  observed  for  divider  separations  of  3.8 
and  8.6  m.  All  the  experiments  reported  here  were  performed  with  divider  separations 
of  3.8  m. 

Figme  40  shows  step-response  measurements  for  the  independent  NBS3  divider  in 
the  horizontal-loop  configuration,  i.e.,  with  no  other  dividers  present  in  the  circuit. 
The  oscillations  aroimd  unit  amplitude  result  from  the  resonance  produced  by  the 
inductance  of  the  high-voltage  lead  and  its  associated  capacitance  along  with  the 
capacitance  of  the  top  ring  of  the  divider.  For  the  NBS3  divider,  this  top  ring 
capacitance  is  smaller,  requiring  a larger  resistance  for  critical  damping  of  the  step- 
response  and  resulting  in  a higher  resonant  frequency  than  that  for  NBSl  which  has 
a larger  top  ring  capacitance.  The  waveforms  in  this  figure  and  in  subsequent  figures 
have  been  displaced  along  the  time  axis  for  better  comparison. 

The  step-responses  of  the  capacitive  dividers  are  shown  in  Figinre  41.  The  pure 
capacitive  divider  had  highly  resonant  behavior  with  a resonant  frequency  of  about 
2.3  MHz  while  the  damped  capacitive  divider  had  a very  slow  rise  time,  but  no  over- 
shoot. Note  that  the  time  scale  of  Figure  41  is  four  times  greater  than  that  of 
Figure  40.  Table  5 summarizes  some  of  the  time  parameters  including  Tn  (experi- 
mental response  time),  T^  (partial  response  time),  and  To  (initial  distortion  time) 
which  were  evaluated  according  to  the  present  standards  [91,  92]  for  the  two  NBS 
resistive  dividers,  NBSl  and  NBS3,  as  a function  of  the  damping  resistor  and  for  the 
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Table  5. 


Rd 

Tn 

Ta 

To 

Overshoot 

(fi) 

(ns) 

(ns) 

(ns) 

(%) 

NBSl 

0 

23 

27 

2.9 

39 

50 

29 

30 

2.4 

28 

150 

33 

36 

2.7 

14 

250 

41 

44 

1.7 

6 

400 

52 

53 

1.8 

2 

NBS3 

0 

19 

17 

2.6 

24 

150 

23 

20 

1.8 

11 

350 

28 

29 

1.3 

3 

550 

36 

36 

0.8 

2 

Damped  Capacitive  Divider 

0 

128 

128 

-0.7 

0 

Pure  Capacitive  Divider 

0 

-23 

46 

0 

49 
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Figure  41.  Step-responses  for  (1)  the  pure  capaeitive  divider,  Rdc  = OQ;  and  (2)  the 
damped  capacitive  divider,  Rddc  = 0 Q. 

damped  and  pure  capacitive  dividers  with  no  damping  resistor. 

Figure  42  shows  examples  of  step-response  measurements  for  selected  configurations 
of  NBSl  with  other  dividers.  For  the  results  shown  in  Fig.  42(a)  no  damping  re- 
sistors were  used  while  in  Fig.  42(b),  NBSl  was  critically  damped.  Independent 
step-response  measurements  are  shown  for  reference.  It  is  clear  that  there  are  in- 
teractions between  the  dividers,  and  that  the  degree  of  interaction  depends  on  the 
damping  resistors.  The  pure  capacitive  divider  had  the  greatest  influence  on  the 
step-response  measured  for  NBSl  for  the  undamped  case.  When  NBSl  was  critically 
damped,  there  was  strong  interaction  observed  with  each  of  the  other  undamped  di- 
viders. If  both  dividers  were  critically  damped,  then  no  interactions  were  observed 
in  the  step  response.  The  distance  between  the  dividers  did  not  greatly  affect  the 
response  until  the  two  dividers  were  located  less  than  1 m apart. 

Undamped  oscillations  in  the  individual  loops  are  the  main  source  of  interaction 
between  the  dividers.  The  divider  system  which  has  the  larger  resonance  amplitude 
has  the  larger  effect  on  the  other  divider.  Among  the  four  dividers,  the  pure  capacitive 
divider  had  the  largest  resonance  amplitude.  The  damped  capacitive  divider  has  no 
resonant  behavior  so  it  was  not  expected  to  and  did  not  greatly  influence  the  step 
responses  of  any  of  the  other  three  dividers. 

Figure  43  shows  examples  of  step-response  measurement  for  NBS3  combined  with 
NBSl  for  different  damping  conditions.  The  step-response  for  NBS3  is  significantly 
distorted  by  an  undamped  NBSl  as  seen  in  Fig.  43(a)  as  is  the  response  for  a critically 
damped  NBS3  (see  Fig.  43(b)). 
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(a) 


Figure  42.  Comparison  of  independent  and  combined  step-responses  for  NBSl. 

(a)  Rdi  = 0 

(1)  Independent  divider 

(2)  Combined  with  pure  capacitive  divider,  R^c  = OQ 

(3)  Combined  with  damped  capacitive  divider,  Rjdc  = Ofl 

(4)  Combined  with  NBS3,  = 0 

(b)  Rdi  = 400n 

(1)  Independent  divider 

(2)  Combined  with  NBS3,  Rds  = 0 

(3)  Combined  with  pure  capacitive  divider,  Rje  = 0f2 

(4)  Combined  with  the  damped  capacitive  divider,  Rjdc  = 0 Q 
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(a) 


Figure  43.  Comparison  of  independent  and  combined  step-responses  for  NBS3. 

(a)  Rd3  = 0 Q 

(1)  Independent  divider 

(2)  Combined  with  NBSl.  Rji  = Ofi 

(3)  Combined  with  NBSl.  Rji  = 400 fi 

(b) R^  = 550Q 

(1)  Independent  divider 

(2)  Combined  with  NBSl.  Rdi  = Ofi 
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Using  damping  resistors  appears  to  be  a simple  and  efficient  way  to  limit  the  inter- 
action caused  by  oscillations,  especially  for  the  divider  measuring  system  which  has 
the  greater  overshoot.  When  both  divider  measuring  systems  had  critical  damping 
resistors,  the  response  from  a combined  measurement  was  nearly  the  same  as  the 
response  observed  from  an  independent  measurement. 

It  may  not  always  be  necessary  to  isolate  the  two  measuring  systems  with  damping 
resistors  because  the  internal  resistance  of  the  impulse  generator  contributes  to  the 
damping  of  oscillations  in  the  measuring  systems.  However,  we  do  need  to  know 
the  interaction  between  two  divider  measuring  systems  when  they  are  simultaneously 
measuring  the  same  signal,  either  a standard  full-  or  chopped-lightning  impulse.  Ob- 
servable interactions  may  introduce  errors  in  a comparative  calibration.  By  studying 
the  step  response  for  the  dividers  connected  for  a simultaneous  measurement,  we  can 
obtain  guidance  as  to  which  configurations  offer  the  smallest  interaction. 


4-2..^  Equivalent  Circuits  of  Measurement  Systems 


Equivalent  circuit  models  were  developed  to  evaluate  the  effects  of  the  interactions 
of  the  dividers  when  an  arbitrary  voltage  waveshape  was  apphed.  The  model  for 
the  independent  divider  NBSl  with  the  step  response  generator  is  illustrated  in  Fig- 
ure 44,  for  which  the  circuit  parameters  were  found  from  either  measurement  or  from 
calculation.  A similar  model  was  developed  for  NBS3.  The  calculated  step  response 
parameters  evaluated  with  these  models  agreed  to  within  5%  of  the  values  evaluated 
from  actual  measurements  for  the  NBSl  divider  while  the  differences  between  the  pa- 
rameters evaluated  from  the  circuit  model  for  the  NBS3  divider  and  from  the  actual 
measurements  were  slightly  larger. 

The  equivalent  circuits  for  the  two  dividers  were  connected  as  shown  in  Figure  45. 
The  inductance  of  the  common  lead  was  estimated  by  calculation  from  the  ringing 
frequency  and  measured  inductance  of  the  horizontal  lead.  A small  coupling  capaci- 
tance between  the  two  dividers  was  also  included  in  the  model.  This  model  was  then 
used  to  simulate  the  simultaneous  measurement  of  the  step  response. 

In  general,  the  calculations  showed  a slightly  larger  interaction  than  those  actually 
observed,  as  seen  in  Figure  46  which  shows  the  measured  step  response  and  that 
computed  from  the  circuit  model.  It  was  shown  through  model  iterations  that  both 
the  inductance  in  the  branch  of  the  connecting  leads  common  to  both  dividers  and  the 
coupling  capacitance  between  them  were  necessary  components  in  the  circuit  model 
to  obtain  agreement  with  the  measured  step  responses,  providing  further  evidence 
that  the  divider  interaction  is  a strong  consideration  in  comparison  measurements. 
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LINE  DIVIDER 


Figure  44.  Equivalent  circuit  model  for  NBSl  divider.  The  circuit  model  includes 
parameters  for  the  4.6  m lead  to  the  relay  and  dc  supply  used  for  the  step  response 
measurements.  The  79  pF  capacitance  is  that  of  the  top  ring  of  the  divider  and  half 
of  the  lead. 


4.2  Interactions  Between  Two  Dividers  Used  in  Comparison  Measurements 


83 


Figure  45.  Equivalent  circuit  used  in  calculation  of  coupled  response  of  NBSl  and 
NBS3.  Blocks  contain  equivalent  circuits  for  NBSl  and  NBS3  (see  Figiue  44).  El- 
ements involved  in  coupling  between  dividers  are  Cc,Lc,  and  Rc-  Rcl  is  a cur- 
rent-limiting resistor. 
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Figure  46.  Comparison  of  calculated  and  measured  NBSl  step  responses  for  simulta- 
neous measurements  using  NBSl  and  NBS3. 

(a)  Rdi  = Rd3  = 0 0 

(1)  Measured 

(2)  Calculated 

(b)  Rdi  =40on,  R^  = on 

(1)  Measmed 

(2)  Calculated 
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TIME,  ps 

Figure  47.  NBSl  divider,  full  standard  lightning  waveform. 

(1)  Input 

(2)  Independent  divider,  R^i  = Ofl 

(3)  Combined  with  NBS3  divider,  Rji  = Rds  = Ofl 

(4)  Combined  with  pure  capacitive  divider,  Rji  = Rdc  = OQ 

(5)  Combined  with  pure  capacitive  divider,  Rji  = 0,  iEjc  = 3500 

(6)  Independent  divider,  Rji  = 3500 

(7)  Combined  with  pure  capacitive  divider,  Rji  = 3500,  = OO 


4.2.5  Standard  Lightning- Waveform  Measurements 


The  interactions  between  dividers  is  clearly  seen  in  the  step  response  measurements 
discussed  in  the  previous  section.  Although  these  interactions  can  be  mitigated 
through  the  use  of  damping  resistors  in  the  connecting  leads  at  the  generator  end, 
it  is  not  always  feasible  to  do  this  in  practice.  Simultaneous  measurements  of  full 
and  chopped  lightning  waveforms  for  diflFerent  combinations  of  dividers  and  circuit 
topologies  were  made  to  observe  the  interactions  and  how  they  can  be  minimized. 

The  output  of  the  NBSl  divider  for  various  combinations  with  other  dividers  and 
damping  resistances  for  a full  standard  lightning  impulse  voltage  of  0.95  ps  rise  time 
and  60  ps  time- to-half- value  are  given  in  Figme  47.  Although  there  are  differences 
in  the  wavefront  measured  with  the  capacitive  divider  when  compared  to  the  inde- 
pendent divider,  the  virtual  front  time  as  evaluated  in  accordance  with  the  stan- 
dards [91,  92]  is  not  significantly  changed:  the  average  for  all  measurements  was 
0.95  ;xs  with  a standard  deviation  of  1%. 
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Measurements  made  with  the  pure  and  damped  capacitive  dividers  for  comparison 
with  NBSl  and  NBS3  are  shown  in  Figure  48.  The  output  of  the  damped  capacitive 
divider  is  significantly  distorted  due  to  this  divider’s  slow  response  time.  The  average 
virtual  front  time,  Ti  found  from  these  measurements  was  1.1  /xs  with  a standard 
deviation  of  0.62%.  This  represents  a 14.4%  front-time  error  which  exceeds  the  10% 
error  limit  for  the  time  parameters  prescribed  by  the  standards  [91,  92].  This  error 
would  only  be  revealed  through  the  comparison  measurements  made  here  and  not 
through  step-response  measurements  since  the  divider’s  measured  unit  step-response 
time  of  129  ns  fell  well  within  the  200  ns  limit  required  by  the  lEC  standard  [91]. 

The  results  for  all  four  dividers  in  independent  and  combined  measurement  of  stan- 
dard lightning  impulses  were  not  significantly  different  except  when  the  pure  ca- 
pacitive divider  was  in  the  circuit.  The  implication  of  this  is  that  the  interactions 
observed  in  the  step-response  measurements  do  not  affect  the  measurement  of  light- 
ning waveforms  except  when  one  of  the  dividers  is  highly  resonant,  such  as  the  pure 
capacitive  divider.  For  standard-lightning  impulses  chopped  at  or  near  the  voltage 
crest,  the  measurements  of  time-to-chop  and  crest  voltages  agreed  very  well  when 
the  NBSl  and  NBS3  dividers  were  simultaneously  used.  The  measurements  for  the 
impulse  voltages  chopped  near  the  crest  did  not  agree  so  well  when  the  capacitive 
divider  was  in  the  circuit.  The  oscillations  present  in  this  case,  as  the  step-response 
measurements  indicate,  give  rise  to  crest  voltage  values  that  are  too  high,  since  the 
time-to-chop  inadvertently  coincided  with  the  crest  of  an  oscillation. 

The  problems  associated  with  evaluating  the  crest  voltage  with  oscillations  present, 
whether  due  to  the  divider  interactions  or  caused  by  the  input  waveform  itself,  are 
clearly  seen  in  Figrue  49(a),  where  measurements  made  with  the  NBSl  divider  in 
combination  with  the  pure  capacitive  divider  are  shown  for  chopped  impulse  voltages. 
The  undamped  pure  capacitive  divider  modifies  the  waveshape  so  that  the  crest 
voltage  is  almost  2%  lower  than  the  individual  measurement  and  the  undershoot  on 
the  voltage  collapse  is  much  larger. 

Figure  49(b)  shows  the  results  from  the  combination  of  the  NBSl  divider  in  the 
critically  damped  configuration  and  the  pure  capacitive  divider.  The  oscillations  are 
nearly  absent  in  this  configuration  and  the  crest  voltages  are  lower  as  a result  of  the 
slower  response  time  of  the  critically  damped  divider. 

For  the  front-chopped  waveforms,  measurements  made  with  the  NBSl  divider  in 
various  combinations  with  other  dividers,  is  shown  in  Figure  50.  The  measurement  of 
the  crest  voltage  depends  upon  the  time-to-chop  for  the  undamped  divider  due  to  the 
severe  oscillations  on  the  wavefront  (Figure  50(a)).  Figure  50(b)  indicates  that  these 
problems  are  not  so  severe  for  the  critically  damped  NBSl  divider.  Only  the  pure 
capacitive  divider  affects  the  measurement  for  this  case  with  the  oscillations  seen  in 
Figure  50(a)  absent,  as  was  anticipated  from  the  step  response  measurements. 


4.2  Interactions  Between  Two  Dividers  Used  in  Comparison  Measurements 


87 


Figure  48.  Response  of  damped  and  pure  capacitive  dividers  to  standard  lightning 
impulses. 

(a)  Damped  capacitive  divider 

(1)  Input  waveform 

(2)  Independent  divider 

(3)  Combined  with  NBSl  divider,  Rji  = 0 

(b)  Pure  capacitive  divider 

(1)  Input  waveform 

(2)  Independent  divider 

(3)  Combined  with  NBSl,  Rji  = OD 

(4)  Combined  with  NBSl,  Rji  = 350 D 
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(a) 


™e^U3 

Figure  49.  NBSl  divider,  chopped  lightning  impulse,  chop  at  crest. 

(a)  Rdi  =01) 

(1)  Input  waveform 

(2)  Independent  divider 

(3)  Combined  with  NBS3  divider,  Rja  = 01) 

(4)  Combined  with  pure  dapacitive  divider,  Rjc  = 01) 

(5)  Combined  with  pure  capacitive  divider,  Rjc  = 5501) 

(b)  Rdi  =4001) 

(1)  Input  waveform 

(2)  Independent  divider 

(3)  Combined  with  NBS3,  Rds  = 01) 


VOLTAGE  VOLTAGE 
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Figure  50.  NBSl  divider,  chopped  lightning  impulse,  chop  at  front. 

(a)  Rdi  = 0 fl 

(1)  Input  waveform 

(2)  Indepenent  divider 

(3)  Combined  with  NBS3  divider,  Rds  = Of! 

(4)  Combined  with  NBS3,  Rja  = 550  fl 

(5)  Combined  with  pure  capacitive  divider,  Rjc  = 0 Q 

(6)  Combined  with  pure  capacitive  divider,  Rjc  = 550 

(b)  R^i  = 400n 

(1)  Input  waveform 

(2)  Independent  divider 

(3)  Combined  with  NBS3,  Rtc  = Ofl 

(4)  Combined  with  NBS3,  R<f3  = 550  f) 

(5)  Combined  with  pure  capacitive  divider,  Rjc  = 0 fl 
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The  strong  interactions  observed  in  the  step-response  measurements  were  less  im- 
portant in  the  measurement  of  standard  full  standard  lightning  impulses  than  in  the 
measurement  of  chopped-lightning  impulses.  This  is  due  to  the  differences  in  the 
frequency  spectra  of  the  two  types  of  waveforms  as  illustrated  in  Figure  51.  The 
spectrum  for  the  full  standard  lightning  waveform  has  a 3dB  bandwidth  of  10  kHz 
and  has  very  small  frequency  components  above  1 MHz,  while  that  for  the  chopped 
lightning  waveforms  has  significant  components  up  to  10  MHz.  The  resonant  oscil- 
lations for  the  divider  systems  had  frequencies  greater  than  2.3  MHz  indicating  that 
for  the  full  standard  lightning  impulse  there  is  not  enough  energy  in  the  frequency 
components  above  1 MHz  to  excite  resonance  even  for  the  undamped  case,  while  for 
the  chopped-lightning  impulse  there  is. 


4.2.6  Conclusions 

The  interaction  between  two  divider  measuring  systems  in  a simultaneous  compar- 
ison measurement  can  produce  distortions  in  the  step-response  waveforms  that  are 
unacceptably  large  according  to  the  relevant  international  standards  [91,  92].  This  is 
especially  of  concern  when  the  resonances  due  to  the  inductances  of  the  connecting 
leads  and  the  stray  capacitances  of  the  leads  and  divider  themselves  are  undamped. 

Although  the  divider  interactions  were  large  in  the  step-response  measurements,  they 
were  smaller  for  the  systems  studied  here  for  full  standard  lightning  impulse  mea- 
surement. The  interactions  produced  errors  in  the  determination  of  the  crest  value, 
time-to-chop,  and  collapse  time  of  chopped-lightning  impulses.  The  dominant  source 
of  error  was  the  lower  frequency  and  higher  amplitude  oscillation  of  the  individual 
undamped  dividers. 

The  differences  in  the  degree  of  interaction  observed  in  the  measurements  of  step- 
response,  full  standard  lightning,  and  chopped-lightning  waveforms  is  explained  quah- 
tatively  through  the  differences  in  their  frequency  spectra  with  respect  to  the  resonant 
frequencies  of  the  dividers.  The  dividers’  resonant  frequencies  are  a few  megahertz 
and  the  normalized  full  standard  lightning  waveform  frequency  spectrum  has  very 
small  components  in  this  range.  The  normalized  chopped-lightning  waveform  has 
larger  amplitude  components  in  this  range  and  has  still  larger  ones  for  the  step- voltage 
waveform.  Thus,  the  interaction  effects  are  most  noticeable  in  the  step-response  mea- 
surements, less  noticeable  in  the  chopped-lightning  impulse  measurements,  and  least 
noticeable  in  the  full  standard  lightning  impulse  measurements. 

The  purpose  of  this  study  has  been  to  observe  the  interactions  of  dividers  used  in 
simultaneous  measurements  in  which  one  divider  is  compared  with  a reference  divider 
to  observe  its  dynamic  response  and  for  evaluation  of  its  scale  factor  (divider  ratio). 
To  perform  this  assessment  acciurately,  it  must  be  demonstrated  that  significant  inter- 
actions do  not  exist.  This  can  be  shown  in  the  following  manner:  first,  by  measuring 
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FREQUENCY,  MHz 


Figure  51.  The  power  spectra  of  (1)  a standard  1.2/60 /zs  lightning  impulse  and  (2)  a 
chopped-lightning  impulse.  The  waveform  was  chopped  at  a voltage  0.9  of  the  crest 
and  the  collapse  time  was  100  ns. 


92 


4.  FAST  TRANSIENT  PHENOMENA 


the  voltage  waveform  of  interest  in  an  undamped  condition.  Then  by  damping  one 
of  the  divider  systems  and  observing  changes  in  measurements  made  with  the  other 
and  by  repeating  the  procedure  with  the  role  of  the  damped  and  undamped  systems 
reversed.  If  significant  interactions  exist,  they  may  be  minimized  through  the  use  of 
damping  resistors  or,  in  some  cases,  through  changing  the  relative  positions  of  the 
divider  systems.  In  any  event,  the  interactions  must  be  proven  neghgible  in  order  to 
acciurately  perform  the  comparison  measurement. 
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